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1 .  Summary 

The  overall  objective  of  this  project  is  to  study  shock  waves  of  the  solar 
terrestrial  medium  in  coronal  and  interplanetary  space.  These  shocks  are 
magnetohydrodynamic  (MHD)  slow  shocks  and  fast  shocks.  Our  primary  interests 
are  to  study  dynamic  problems  related  to  shocks  associated  with  coronal  mass 
ejections  (CMEs),  including  the  formation  of  shocks  in  coronal  environment, 
the  transition  of  slow  shocks  to  fast  shocks,  their  evolution  and  their 
interaction  with  the  solar  wind  in  coronal  and  interplanetary  space. 

During  the  past  three  years  of  this  grant,  we  have  produced  the  following 
publications  : 


[Paper  1]  Whang,  Y.  C.,  Slow  Shocks  and  Their  Transition  to  Fast  Shocks  in  the 
Inner  Solar  Wind,  J.  Geophys,  Res.,  92,  4349-4356  (1987). 

[Paper  2]  Whang,  Y.  C.,  Evolution  of  Interplanetary  Slow  Shocks,  J.  Geophys. 
Res.,  93,  251  -255  (1988).  '  '  _ 

[Paper  3]  Whang,  Y.  C.,  Forward-Reverse  Shock  Pairs  Associated  with  Coronal 
Mass  Ejections,  J.  Geophys.  Res.,  93 ,  5897-5902  (1988). 

[Paper  4]  Whang,  Y.  C.,  CME  Associated  Forward-Reverse  Shock  Pairs,  Adv, 

Space  Tes.,  9,  (4)81-86,  (1989). 

[Paper  5]  Waang,  Y.  C.,  Formation  of  Slow  Shock  Pairs  Associated  With  Coronal 
Mass  Ejections,  in  press,  in  the  Proceedings  of  the  Chapman  Conference  on 
the  Physics  of  Magnetic  Flux  Ropes,  ACU  Publication  (1989). 

In  the  following  sections,  we  summarize  some  highlights  of  these  studies. 


2 .  Formation  of  Slow  Shock  Pairs  Associated  With  CMEs 

In  order  to  understand  the  basic  physics  involved  in  the  shock  formation 
process,  we  carry  out  a  series  of  MHD  simulations  to  show  that  MHD  shocks 
formed  in  low-beta  plasmas  in  the  coronal  environment  are  slow  shocks.  The 
simulation  model  assumes  that  flow  properties  are  functions  of  x  and  t  only. 
The  x-direction  is  normal  to  the  wave  front  in  the  forward  direction  of  the 
solar  wind  flow  and  the  magnetic  field  vector  is  parallel  to  the  xy  coordinate 
plane.  The  time-dependent  MHD  flow  governed  by  a  system  of  hyperbolic, 
nonlinear  equations  is  studied  using  the  method  of  characteristics.  Shocks  are 
treated  as  surfaces  of  zero  thickness.  The  Rankine -Hugoniot  solution  is  used 
to  calculate  the  jumps  in  flow  properties  at  all  shock  crossings. 

The  numerical  solution  shown  in  Figure  1  demonstrates  the  formation  of  a 
forward- reverse  slow  shock  pair  in  a  coronal  environment.  The  initial 
condition  represents  the  impact  of  a  high-speed  (600  km/s)  mass  ejecta  on  a 
low-speed  (100  km/s)  ambient  solar  wind.  At  t  -  0  the  flow  velocity  changes 

mass  ejecta) 
tream  (ambient 
.  The  flow  is 

field  aligned  in  the  rest  frame  of  reference,  0  -  0.1,  Bx  -  0.8  Gauss,  and  6  - 


smoothly  over  a  distance  of  about  200,000  km.  On  the  upstream  ( 
side  -  600  km/s  and  ni  -  8.8  x  10^  protons/cm^.  On  the  downs 
solar  wind)  side,  U  -  100  km/s  and  n  -  6.0  x  10^  protons/cm 
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The  simulation  demonstrates  that  the  difference  in  momentum  flux  is  the 
driving  force  for  the  formation  of  shock  pairs.  The  momentum  impact  compresses 
the  plasma  near  the  front  of  the  CME.  Large  pressure  disturbances  propagate  in 
both  the  forward  direction  relative  to  the  ambient  solar  wind  and  the  reverse 
direction  relative  to  the  ejecta  flow.  The  pressure  fronts  steepen  to  from 
forward  and  reverse  MHD  shocks.  The  flow  speed  profiles  calculated  at  a  time 
interval  of  30  seconds  (Figure  1)  show  the  gradual  formation  of  the  two  shock 
fronts.  The  resulting  shocks  consist  of  a  forward  slow  shock  and  a  reverse 
slow  shock.  No  fast  shocks.  The  shock  pair  takes  about  10  minutes  to  grow  into 
a  fully  developed  state. 

The  disturbances  generated  by  the  momentum  impact  also  disturb  the  flow 
field  outside  the  interaction  region  bounded  by  the  forward- reverse  slow  shock 
pair.  But  their  influence  is  limited  to  the  region  bounded  by  a  forward  fast 
mode  characteristic  on  the  downstream  side  and  by  a  reverse  fast  mode 
characteristic^  or  the  am  side.  The  amplitudes  of  disturbances  outside 

the  slow  shocks  are  negligibly  small  compared  with  jumps  in  flow  properties 
across  the  slow  shocks. 


U,  =  600  km/s  U0  =  1 00  krn/s  N,/N0  =  1.I46  0=0.10=15° 
-11  3  5  7  9 


Fig.  1.  The  flow  speed  profiles  calculated  at  a  time  interval  of  30  seconds 
showing  the  gradual  formation  of  the  two  shock  fronts. 
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3 .  A  Parametric  Survey  of  MHD  Slow  Shocks 


We  have  an  exact  analytical  method  to  calculate  the  solutions  of  MHD 
shocks  as  functions  of  three  dimensionless  upstream  parameters:  uhe  shock 
Alfven  number  A  -  Un/(a  cos  8),  the  shock  angle  8,  and  the  plasma  p  value. 
(Here  Un  is  the  normal  component  of  the  relative  shock  speed,  a  the  Alfven 
speed,  8  the  acute  angle  between  the  shock  normal  and  the  magnetic  field,  and 
P  the  ratio  of  the  thermal  pressure  p  to  the  magnetic  pressure  B^/8 n.)  A  shock 
is  a  slow  shock  if  the  shock  Alfven  number  A  is  less  than  1  and  is  a  fast 
shock  if  A  is  greater  than  1.  The  magnetic  field  and  the  shock  angle  decrease 
across  a  slow  shock,  they  increase  across  a  fast  shock. 

The  parametric  survey  of  slow  shock  solutions  are  presented  in  a  three 
dimensional  A,0  ,/J-space.  The  domain  of  slow  shock  solutions  in  the  A  ,5- 
parameter  space  shrinks  as  the  plasma  P  value  decreases.  A  drastic  change  in 
the  domain  of  solution  takes  place  when  P  is  between  0.1  and  1.  The  change 
becomes  very  small  at  very  large  or  very  small  p.  When  the  shock  Alfven  number 
approaches  one,  slow  shocks  exist  at  all  shock  angles  and  at  all  p  values.  At 
any  given  value  of  P,  a  whole  range  of  changes  in  B2/B1,  from  0  to  1,  may  take 
place.  On  the  other  hand  the  range  of  change  in  thermodynamic  properties 
strongly  depends  on  the  p  value.  In  Figure  2  of  Paper  2  shows  the  density 
ratio,  P2/P\<  i-n  t*ie  three  dimensional  A, 6, ^-parameter  space.  (A  similar  plot 
has  been  made  for  the  pressure  ratio.)  Slow  shocks  covering  a  wide  range  cf 
density  jumps  exist  in  the  low  p  region  (where  p  <  0.1).  On  the  other  hand, 
only  slow  shocks  with  weak  density  jumps  exist  in  the  high  p  region  (where  p  > 
1).  Therefore,  in  the  coronal  space  where  the  p  value  of  the  solar  wind  is  in 
the  order  of  0.1,  all  physical  properties  may  jump  across  a  slow  shock  over  a 
wide  range  of  magnitudes.  Near  1  AU  the  plasma  p  value  is  in  the  order  of  1, 
the  jumps  in  thermodynamic  properties  across  a  slow  shock  must  be  small  but 
the  change  in  the  magnetic  field  is  not  necessarily  small. 


4.  Transition  of  Slow  Shocks  to  Fast  Shocks 


Whang  [Paper  1]  first  studied  the  possible  existence  of  forward  slow  shocks 
preceding  some  CMEs  and  their  global  geometry.  When  the  speed  of  ejected  CME 
material  relative  to  the  ambient  solar  wind  exceeds  the  local  magnetoacoustic 
speed,  an  MHD  shock  front  forms  at  the  leading  edge  of  the  compressed  ambient 
plasma  shell.  In  coronal  space,  the  Alfven  speed  is  of  the  order  of  1000  km/s. 
If  a  forward  MHD  shock  forms  preceding  the  compressed  ambient  plasma  shell 
associated  with  the  CME,  the  shock  Alfven  number  can  be  less  than  1  and  the 
shock  can  be  a  forward  slow  shock.  The  large-scale  global  surface  of  a  forward 
slow  shock  should  have  a  bow-shaped  surface  with  its  nose  facing  the  sun 
because  the  slow-mode  MHD  waves  and  MHD  shocks  are  upstream-inclined  [Sears 
and  Resler,  1961],  The  relative  shock  speed  varies  from  point  to  point  on  the 
curved  surface  of  an  MHD  shock.  The  slow  shock  is  relatively  stronger  at  the 
nose.  The  shock  becomes  weaker  on  the  flank  and  asymptotically  approaches 
slow-mode  MHD  waves. 

Sheeley  et  al.  [1985]  have  demonstrated  a  close  association  between 
interplanetary  shocks  and  CMEs.  They  reported  that  with  very  few  exceptions 
interplanetary  shocks  detected  at  Hellos  A  are  produced  by  CMEs  observed  from 
Solwind  coronagraph.  Richter  et  al.  [1985]  reported  a  very  certain 
identification  of  a  forward  slow  shock  at  0.31  AU.  Hundhausen  et  al.  [1987] 
studied  CMEs  observed  from  the  coronagraphs  aboard  Skylab  and  on  the  Solar 
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Maximum  Mission  satellite.  The  observed  flattening  of  the  tops  of  some  mass 
ejection  loops  suggests  the  possible  existence  of  traveling  forward  slow 
shocks  immediately  in  front  of  some  CMEs.  The  observed  speeds  for  the  outward 
motion  of  CMEs  show  that  the  speed  at  which  a  CME  is  overtaking  the  ambient 
coronal  plasma  is  less  than  the  Alfven  speed  but  greater  than  the  sound  speed 
for  many  (probably  most)  observed  CMEs.  This  supports  the  possible  existence 
of  forward  slow  shocks. 

The  Alfven  speed  of  the  interplanetary  plasma  decreases  at  increasing 
heliocentric  distance  in  the  inner  solar  wind.  As  a  forward  slow  shock  travels 
outwards  from  the  sun,  the  decrease  in  the  Alfven  speed  causes  the  shock 
Alfven  number  of  a  forward  slow  shock  to  become  greater  than  1  at  increasing 
heliocentric  distances  and  the  shock  evolves  from  a  slow  shock  into  a  fast 
shock.  The  large-scale  global  surface  of  a  traveling  forward  fast  shock  has  a 
bow-shaped  surface  with  its  nose  facing  away  from  the  sun. 

Whang  [Paper  1]  has  a  theory  to  explain  the  transition  of  traveling  »iov 
shocks  into  traveling  fast  shocks  as  shown  in  Figure  6  of  Paper  1.  The  Alfven 
speed  of  the  interplanetary  plasma  decreases  at  increasing  heliocentric 
distance  in  the  inner  solar  wind.  As  the  forward  slow  shock  moves  outwards 
from  the  sun,  the  decrease  in  the  Alfven  speed  causes  the  shock  Alfven  number 
to  reach  1  near  the  nose  of  the  shock  surface,  the  transition  from  a  forward 
slow  shock  to  a  forward  fast  shock  begins  to  take  place  at  a  point  of  the 
shock  surface  where  the  shock  angle  equals  0°.  The  onset  of  transition  must 
occur  in  the  region  of  interplanetary  space  where  ,5  <  2/y  or  1.2.  At  the  onset 
of  the  transition  the  forward  slow  shock  smoothly  converts  to  a  system 
consisting  of  a  forward  slow  shock,  a  very  weak  forward  rotational 

discontinuity  (an  Alfven  wave),  and  a  very  weak  forward  fast  shock  (a  fast  MHD 

wave) . 

During  the  transition,  the  system  consists  of  a  slow  shock,  a  fast  shock, 
and  a  rotational  discontinuity.  The  three  surfaces  of  discontinuity  intersect 
along  a  closed  transition  line.  As  the  system  moves  outward  from  the  sun  the 
transition  line  moves  laterally  across  the  field  lines,  the  area  enclosed  by 
the  closed  loop  of  the  transition  line  expands  continuously,  the  fast  shock 
grows  stronger,  and  the  slow  shock  becomes  weaker.  Eventually,  the  slow  shock 
diminishes  and  the  entire  system  evolves  into  a  forward  fast  shock.  Such  a 
system  may  have  been  observed  from  ISEE  1  and  3.  Kennel  et  al.  [1984a,  b] 

reported  that  on  November  12,  1978,  ISEE  1  and  3  observed  a  forward  fast  shock 

associated  with  a  pair  of  flares  which  occurred  some  48  hours  earlier.  The 
shock  was  followed  by  a  strong  rotational  discontinuity  at  approximately  40  Rg 
downstream  of  the  shock.  This  sequence  of  events  resembles  the  forward 
transition  system. 


5.  Evolution  of  CME  Associated  Shock  Pairs 


Recently,  Whang  [Paper  3]  suggests  that  for  some  CMEs  their  interaction 
with  the  ambient  solar  wind  can  produce  a  forward-reverse  shock  pair  as  shown 
in  Figure  1  of  Paper  3.  We  may  call  the  ejected  material  the  CME  plasma  which 
is  separated  from  the  solar  wind  plasma  in  the  background  corona  by  an 
interface.  The  large  momentum  of  the  high-speed  mass  ejecta  exerts  an 
excessive  pressure  on  the  ambient  plasma  shell  in  a  narrow  region  in  front  of 
the  CME.  The  high-speed  mass  ejecta  causes  a  sudden  increase  in  the  pressure 
of  the  plasma  near  the  top  of  the  CME  on  both  sides  of  the  interface. 
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In  response  to  the  large  pressure  in  the  compressed  ambient  plasma  shell, 
the  CME  plasma  immediately  behind  the  interface  is  decelerated,  compressed  and 
heated.  This  compressed  CMF,  plasma  is  confined  in  a  narrow  region  bounded  by 
the  interface  on  the  forward  side  and  by  a  pressure  front  on  the  reverse  side. 
The  front  of  the  compressed  CME  plasma  propagates  in  the  reverse  direction 
relative  to  the  ejecta  flow,  it  may  steepen  to  from  a  reverse  shock. 

The  ejected  CME  plasma  carrying  a  strong  magnetic  field  is  a  low-/9  plasma 
in  which  the  slow  mode  magnetoacoustic  speed  Cg  is  much  less  than  the  Alfven 
speed  a.  The  reverse  shock  associated  with  CME  can  be  a  slow  shock  if  the 
normal  component  of  the  relative  shock  speed  Un  is  greater  than  Cs  but  less 
than  a  cos  6.  Whang  estimated  that  the  reverse  shock  can  be  a  slow  shock 
whether  the  forward  one  is  a  slow-mode  or  a  fast- mode  shock.  As  a  slow  shock 
pair  travels  outward  from  the  Sun,  the  decrease  in  the  Alfven  speed  with 
increasing  heliocentric  distance  causes  the  slow  shocks  to  evolve  into  fast 
shocks . 

MacQueen  [1980]  suggested  that  CMEs  do  not  carry  out  distended  field  very 
far  to  interplanetary  space,  but  rather  are  subject  to  a  reconnection  process. 
He  proposed  that  CME  loops  must  magnetically  disconnect  from  the  sun  and  form 
a  closed  magnetic  structure. 

Burlaga  et  al.  [1981]  analyzed  the  magnetic  field  configuration  behind 
three  shocks.  They  found  a  systematic  configuration  of  the  magnetic  field  and 
called  it  the  magnetic  cloud.  Klein  and  Burlaga  [1982]  identified  45  clouds  in 
interplanetary  data  obtained  near  1  AU  between  1967  and  1978.  They  compared 
the  physical  properties  of  magnetic  clouds  with  those  of  CMEs  and  suggested  an 
association  of  some  clouds  with  disconnected  magnetic  structures  of  CMEs.  They 
identified  a  magnetic  cloud  as  an  interplanetary  structure  in  which  the 
magnetic  field  strength  is  higher  than  average,  the  magnetic  fi^ld  direction 
rotates  monotonically  through  a  large  angle,  the  temperature  is  low,  and  the 
plasma  P  is  significantly  lower  than  1.  A  good  case  of  an  association  between 
a  CME  observed  by  Solwind  and  a  magnetic  cloud  observed  at  0.54  AU  by  Helios  1 
two  days  later  was  presented  by  Burlaga  et  al.  [1982], 

In  summary,  two  important  evolutions  of  the  shock  pairs  associated  with 
CMEs  may  occur  in  interplanetary  space.  First,  a  CME  loop  eventually 
disconnect  from  the  sun  to  form  a  closed  magnetic  structure.  The  disconnected 
bubble  manifests  as  a  magnetic  cloud  in  interplanetary  space.  Second,  as  the 
CME  associated  forward-reverse  shock  pair  moves  outwards  in  interplanetary 
space,  they  evolve  into  a  pair  of  fast  shocks.  The  reverse  shock  first 
propagates  within  the  magnetic  cloud,  but  eventually  the  reverse  shock  must 
exit  the  cloud  to  became  detached  from  the  magnetic  cloud.  Then  both  the 
forward  and  the  reverse  shocks  propagate  in  the  ambient  solar  wind.  Therefore, 
the  interplanetary  consequence  of  some  CMEs  should  consist  of  a  magnetic  cloud 
accompanied  by  a  fast  shock  pair.  The  forward  fast  shock  propagates  preceding 
the  cloud  and  the  reverse  fast  shock  is  either  within  or  behind  the  cloud. 

Klein  and  Burlaga  [1982]  reported  that  about  one  third  of  clouds  identified 
near  1  AU  between  1967  and  1978  are  preceded  by  a  shock.  There  are  several 
examples  involving  the  observation  of  a  shock  pair  associated  with  a  CME  or 
magnetic  cloud.  Lepping  et  al.  [1988]  reported  that  a  large  magnetic  cloud 
accompanied  by  a  shock  pair  was  observed  on  October  31  and  November  1  of  1972 
from  IMP  7.  The  event  was  associated  with  a  unusually  longlasting  solai.  flaie 


which  occurred  some  49  hours  earlier.  Gosling  et  al.  [1988]  have  identified 
two  CME  associated  shock  pairs  at  1  AU  from  ISEE  3.  Another  event  involving  a 
shock  pair  and  a  magnetic  cloud  was  observed  on  August,  1982,  from  Voyager  2 
at  10.3  AU  [Burlaga  et  al.,  1985].  This  observation  shows  that  for  a  long- 
lived  cloud,  the  reverse  shock  had  enough  time  to  propagate  through  the  cloud, 
exited  the  rear  of  the  cloud,  and  appeared  closely  behind  the  cloud. 
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Abstract.  The  jump  conditions  of  MHD  shocks 
may  be  directly  calculated  as  functions  of  three 
upstream  conditions:  the  shock  Alfven  number 
based  on  the  normal  component  of  the  relative 
shock  speed,  the  shock  angle,  and  the  plasma  $ 
value.  The  shock  Alfven  number  is  less  than  1  for 
a  slow  shock  and  greater  than  1  for  a  fast  shock. 
A  traveling,  forward  shock  can  be  a  slow  shock  in 
coronal  space,  where  the  Alfven  speed  is  of  the 
order  of  1000  km/s.  The  surface  of  a  forward  slow 
shock  has  a  bow-shaped  geometry  with  its  nose 
facing  toward  the  sun.  The  decrease  in  the  Alfven 
speed  at  increasing  heliocentric  distance  causes 
the  shock  Alfven  number  of  a  forward  slow  shock 
to  become  greater  than  1,  and  the  shock 
eventually  evolves  from  a  slow  shock  into  a  fast 
shock.  During  the  transition  the  shock  system 
consists  of  a  slow  shock,  a  fast  shock,  and  a 
rotational  discontinuity.  They  intersect  along  a 
closed  transition  line.  As  the  system  moves 
outward  from  the  sur,  the  area  enclosed  by  the 
transition  line  expands,  the  fast  shock  grows 
stronger,  and  the  slow  shock  becomes  weaker. 
Eventually,  the  slow  shock  diminishes,  and  the 
entire  shock  system  evolves  into  a  forward  fast 
shock . 

1.  Introduction 

Interplanetary  traveling  shocks  are  associated 
with  c  .’ona  mass  election:,  and  flares  IDrver 
1975;  Hildner,  19/7;  MacQueen,  1980;  Sheelev  et 
al..  1985],  If  the  speed  of  the  ejected  material 
relative  to  the  ambient  solar  wind  exceeds  the 
local  magnetoacoustic  speed,  a 
magnetohydrodynamic  (MHD)  shock  front  forms  at 
the  leading  edge  of  the  compressed  ambient  plasma 
shell  Th®  shock  process  plays  an  important 

role  in  influencing  solar  terrestrial  teitilicr.o. 

The  solutions  of  MHD  shocks  may  be  expressed 
as  functions  of  three  upstream  conditions:  the 
shock  Alfven  number  based  on  the  relative  shock 
speed  A  -  Un/(a  cos  9),  the  shock  angle  6,  and 
the  plasma  value.  (Here  U  is  the  normal 
component  of  the  relative  shock  speed,  a  the 
Alfven  speed,  9  the  acute  angle  between  the  shock 
normal  and  the  magnetic  field,  and  8  the  rat.  lo  of 
the  thermal  pressure  p  to  the  magnetic  pressure 
B  /8tt.)  Two  kinds  of  MHD  shocks,  the  fas*  (<ud  the 
slow  shocks,  may  exist.  A  shock  is  a  slow  shock 
if  the  shock  Alfven  number  A  is  less  than  1  and 
is  a  fast  shock  if  A  is  greater  than  1  . 

A  traveling,  forward  MHD  shock  can  he  a  slow 
shock  in  coronal  space  where  the  Alfven  speed  i s 
of  the  order  of  1000  km/s.  The  Alfven  speed 
decreases  at  increasing  heliocentric  distance  in 
the  inner  solar  wind.  The  decrease  in  the  Alfven 
speed  causes  the  shock  Alfven  number  of  a  forward 
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slow  shock  to  become  greater  than  1  at  i  nr  re.*s :  r.g 
heliocent  ric  distances,  and  the  shock  «-vrr:'  !v 
evolves  trom  a  slow  shock  into  a  fast  .chock.  A 
theoretical  prediction  for  the  global  geo met rv  of 
traveling  slow  shocks  and  ‘.heir  transition  to 
fast  shocks  was  first  presented  at  Spring 

Meeting  of  the  American  Geophysical  Union  [Wh.ing. 
1986a],  Recently,  we  have  carried  out  a  4iuTt-ric.il 
solution  to  demonstrate  the  transition  process 
This  paper  contains  a  formal  presentation,  of 
these  researches. 

Theoretical  investigations  of  MHD  slow  shocks 
in  space  plasma  have  been  reported  by  Neubauer 
[1976],  Rosenau  and  Suess  [19/7],  Whang  [19S2, 
1986a,  b  ]  ,  Hada  and  Kennel  [1985].  and  Edmistor. 
and  Kennel  [1986].  Very  few  of  the  observed 
interplanetary  shocks  are  slow  shocks. 
Observations  of  interplanetary  slow  shocks  near  1 
AU  have  been  reported  by  Chao  and  Olbert.  [19’0[ 
and  by  Burlaga  and  Chao  [1971].  However, 
observations  of  slow  shocks  by  Helios  at  a 
heliocentric  distance  of  0.31  AU  [Richter  et  al., 

1985]  and  by  ISEE  3  beyond  100  earth's  radii  in 
the  geomagnetic  tail  [Feldman  et  al.,  1985] 
greatly  increase  our  interest  in  slow  shocks.  A 
better  understanding  of  the  MHD  slow  shock 
process  may  have  direct  applications  to  the 
magnetic  field  reconnection  process  and  to  shock 
processes  in  coronal  space,  in  interplanetary 
space,  and  in  geomagnetic  tails. 

Recently,  A.  J.  Hundhausen  et  al.  (Do  slow 
shocks  precede  some  coronal  mass  ejection?, 
su’  m  i  1 1  ed  to  Juu;  na  1  of  Geophvs  i  ca  1  Re  sea  ren, 

1986)  (hereafter  referred  to  as  HHL)  studied  the 

outward  motion  for  coronal  mass  ejections 
observed  by  the  Skylab  eoronagraph  and  the  SMM 
coronagraph.  The  measured  speeds  of  coronal  mass 
ejections  indicate  that  many  mass  ejections  move 
in^  rhr  v-  cictrtr  than  the 

gasdynamic  sound  speed  hut  slower  than  the  Alfven 
speed.  This  has  been  interpreted  as  a  support  for 
the  possible  formation  of  slow  MHD  shocks 
immediately  in  front  of  some  coronal  mass 

e  j  ec  t i ons . 

2.  Interplanetary  Slow  Shock  and  Its  Trans  on¬ 
to  Fast  Shock 

Two  kinds  of  MHD  shocks,  the  last  and  the  s  low- 
shocks.  may  exist:  in  coronal  and  interplanetary 
space.  The  plasma  density,  the  thermal  pressure, 
and  the  total  pressure  (the  sum  of  the  thermal 
and  the  magnetic  pressures)  all  increase  across 
an  MHD  shock.  The  magnetic  field  and  the  shock 
angle  decrease  across  a  slow  shock  and  increase 
across  a  fast  shock  The  density  ratio  ei.-Vj  is 
often  used  as  a  measure  of  the  shock  strength. 
Here  the  subscript  1  denotes  th*  flow  condit  i on 
upstream  of  a  shock,  and  ?  denotes  the  flow 
condit  ion  downstream.  Fipures  1  and  ?  show  tin* 
constant  contours  foi  tne  density  ratio  and  the 
magnetic  field  ratio  for  slow  and  fa s ’  shocks, 
respectively  as  functions  of  three  dimensionless 
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Alfven  Number  =  Un/(a  cos  9) 

Fig.  1.  The  constant  contours  for  the  density 
ratio  and  the  magnetic  field  ratio  of  slow 
shocks.  The  jump  conditions  of  MHD  shocks  can  be 
directly  calculated  as  functions  of  three 
dimensionless  upstream  parameters:  the  shock 
Alfven  number  A  -  Un^  /{a^  cos  tf^),  thr  shock 
angle  9 ^ ,  and  ft-^.  The  magnitude  of  the  magnetic 
field  decreases  across  a  slow  shock. 


upstream  parameters:  the  shock  Alfven  number  A  - 
Unl  /(a:  cos  sp,  Che  shock  angle  and  fiv 

In  coronal  space  the  Alfven  speed  is  of  the 
order  of  100O  k m/s;  an  MHD  shock  in  the  inner 
solar  wind  is  a  slow  shock  if  the  shock  Alfven 
number  A  is  less  than  1.  An  interesting  feature 
of  interplanetary  slow  shocks  is  that  the  surface 
of  a  forward  slow  shock  should  have  a  bow-shaped 
surface  with  its  nose  facing  the  sun.  The  sketch 
on  the  left  panel  of  Figure  3  shows  the  predicted 
geometry  of  a  forward  slow  shock  on  the  helio¬ 
graphic  equator.  The  solar  wind  plasma  enters  a 
shock  from  its  front  side.  The  front  side  (the 
upstream  side)  of  a  forward  shock  is  the  side 
farther  away  from  the  sun  and  back  side  (the 
downstream  side)  is  the  side  closer  to  the  sun. 

An  important  ch irac terist ic  of  slew  sh  (and 

slow  MHD  waves)  is  that  they  are  upstream  - 
inclined  [Sears  and  Resler,  1961].  Therefore, 
locking  at  the  forward  slow  shock  from  the  sun, 
one  sees  the  convex  back  side  surface  of  the  slow 
shock.  The  relative  shock  speed  varies  from  po{nt 
to  point  on  the  curved  surface  of  an  MHD  shock. 
The  slow  shock  is  relatively  stronger  at  the 
nose.  The  shock  becomes  weaker  on  the  flank  and 
asymptotically  approaches  to  slow  mode  MHD  waves. 
The  shock  angle  6  decreases  across  the  slow  shock 
by  a  deflection  angle  6.  HHL  pointed  out  that 
the  upstream  inclined  nature  of  slow  MHD  shocks 
is  relative  to  the  magnetic  field  geometry;  thus 
in  a  diverging  radial  background  field,  a  large- 
scale  traveling  slow  shock  may  have  a  rather  flat 
configuration.  They  used  the  observed  flattening 
of  the  tops  of  some  mass  ejection  loops  as  they 
move  outward  through  the  corona  to  support  the 
possible  existence  or  traveling  •*»’.:  shocks. 

In  the  inner  heliosphere  the  Alfven  speed 
decreases  at  increasing  heliocentric  distance.  As 
a  forward  slow  shock  travels  outward  from  the 
sun.  the  decrease  in  the  Alfven  speed  causes  the 
shock  Alfven  number  of  a  forward  slow  shock  to 
reach  1  near  the  nose  of  the  shock  surface,  and 
the  transition  from  a  slow  shock  to  a  fast  shock 


begins  to  take  place.  The  shock  angle  9  -  0  at 
the  initial  transition  point.  The  onset  of 
transition  can  occur  in  the  region  where  ft  e  2/* y 
oi  i .  ? .  Here  7  is  the  ratio  of  speritic  heats. 

All  calculations  carried  out  in  this  paper  use  7 

-  5/3  for  a  fully  ionized  plasma  Note  that  when 

3  -  2/7,  the  Alfven  speed  equals  the  gasdynamic 
sound  speed  When  ft  <  1.2,  the  phase  speed  of  the 
fast  MHD  wave  equals  the  Alfven  speed  at  6  -  0 
The  plots  ir.  Figure  2  denicr.stiate  that  near  A 

-  1,  fast  MHD  shocks  exist  when  ft  <  1.2  and  but 
not  when  ft  >  1.2 

At  greater  heliocentric  distances  the  forward 
slow  shock  may  eventually  evolve  into  a  forward 
fast  shock.  Looking  at  a  forward  fast  shock  from 
the  sun,  the  back  side  surface  of  the  forward 
traveling  fast  shock  is  concave  because  the  fast 
mode  MHD  waves  and  MHD  shocks  are  downstream- 
inclined.  Thus  the  large-scale  global  surface  of 
a  traveling  forward  fast  shock  has  a  bow-shaped 
surface  its  nose  facing  away  from  the  sun  as 

shown  in  the  right  panel  of  Figure  3.  The 
deflection  angle  of  the  magnetic  field  represents 
the  increase  of  the  shock  angle  across  a  fast 
shock . 

During  the  transition  the  shock  system 
consists  of  a  slow  shock,  a  fast  shock,  and  a 
rotational  discontinuity  which  may  be  referred  to 
as  a  slow-fast  shock  system.  We  will  continue  to 
discuss  the  slow-fast  shock  system. 

3.  Jump  Conditions  of  MHD  Shocks 

We  obtain  a  direct  and  simple  solution  for  the 
jump  conditions  of  MHD  oblique  shocks.  The  jump 
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Alfven  Number  =  Un/(a  cos  9) 


Fig.  2.  The  constant,  contours  for  the  density 
ratio  and  the  magnetic  field  ratio  of  fast 
shocks.  The  magnitude  of  the  magnetic  field 
increases  across  a  fast  shock  Near  A  -  1,  fast 
shocks  exist  when  ft  ^  <  1.2  but  not  when  ft ^  >  l.  2. 
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Fig  3.  The  large-scale  global  surface  of  a 
traveling  forward  slow  shock  should  have  a  bow  • 
shaped  surface  with  its  nose  facing  the  sun  as 
shown  on  the  left,  panel.  The  shock  angle  6 
decreases  across  the  slow  shock  bv  a  deflection 
angle  6.  As  a  forward  slow  shock  travels  outward 
from  the  sun.  the  decrease  In  the  Alfven  speed 
causes  the  shock  Alfven  number  of  a  forward  slow 
shock  to  reach  1  near  the  nose  of  the  shock 
surface,  and  the  transition  from  a  slow  shock  to 
a  fast  shock  begins  to  take  place.  At  large 
heliocentric  distances  the  slow  shovk  is  totally 
evolved  into  a  fast  shock  The  large-scale  global 
surface  of  a  traveling  forward  fast  shock  has  a 
bow-shaped  surface  with  its  nose  facing  away  from 
the  sun.  as  shown  on  the  right  panel.  The  shock 
angle  increases  across  a  fast  shock 


-  s  l,  > .  Bt,rf  -  rr.  b.  > '  nt  ; 

and 

Rn'V  -  rnB-- 

From  equations  (4)  and  ■  i ,  we  can  obtain 

£*xp\}  U  B  '  -  B  B  ’ 
h  n 1  n  r  •  n  t 

Making  use  of  equations  ( *> )  ,  '  '  ,  and  i'  8 )  .  we 

w.-ite  the  energy  equation  (b)  as 


yp/r):,y  -  1  - ' 


We  shall  show  that  for  given  upstream 
conditions,  from  i  ? )  .  < .3  )  .  (  8  >  .  and  i 4  >  w e  c  a  n 
obtain  the  exact  solution  of  the  mum  shook 
equations  for  p-s.  B^-,.  p  ,  and  U_-,  as  functions 
of  three  dimensionless  upstream  parameters  the 
shock  Alfven  number  A,  the  shock  angle  9  ^  ,  and 
,;h  Note  that  these  solutions  are  independent  o 
the  velocity  components  l*t  and  t's-  [l\'  can  he 
direct  lv  calnilated  from  14)  The  shock  .solvit  in 
remain  unchanged  if  an  arbitrary  l’s  is  added  to 
bo t  h  sides  of  r he  s hoc k 

It  is  convenient  to  denote  the  density  ratio 
bv 


V  -  fx/V;  -  VniA'n2 

Th.*-n  from  (?)  and  i  /  )  we  obtain 

n  -  aV,/!(a2  -  UBtl  ♦  B.,! 


conditions  of  MHD  shocks  can  he  calculated  in 
terms  of  some  shock  parameters,  which  are 
sometimes  defined  in  terms  of  physical  properties 
on  both  sides  of  the  shock  This  method  of 
solution  has  all  shock  parameters  defined  in 
terms  of  the  upstream  conditions  only. 

It  is  convenient  to  formulate  the  shock 
solution  using  an  n,  t . s - coordi nate  system 
attached  to  the  shock  surface.  Let  n  be  the  unit 
vector  normal  to  the  shock  surface  and  pointing 
in  the  direction  of  mass  flow  Since  the  shock 
normal  n  and  the  magnetic  field  vectors  on  both 
sides  of  a  shock.  B|  and  B  , ,  are  eoplanar,  this 
coplanarity  plane  is  chosen  as  the  nt  -coordinate 
plane.  This  means  that  Bs  -  0  on  both  sides  of 
the  s  hoc  k . 

We  use  a  pair  of  brackets  to  denote  the  jump 
of  a  physical  quant itv  Q  across  the  shock,  e  g  . 

f q ;  -  0,  - 

Let  U  be  the  flow  velocity  in  a  frame  of 
reference  attached  to  the  shock  surface.  At  anv 
instant  the  jumps  in  plasma  and  magnetic  fields 
across  an  MMD  shuck  must  satisfy  the  following 
jump  equation*;  'Heifer, 

'  B  ’  =  0  .it 

n  ' 

1  p\"  '  -  O  (  ?  i 

n 

c  \ '  '  l  i  ’  *-  ;  p  :  ♦  M  /  8  IT  ')  B *  :  -  0  (  3  > 

n  •  n j  1  *  1  t 

;'cl'.rl:t  i  -  RnB:  *  1  4  1 

U's  -  "  C-O 

-  1)  i 


and 

V  -  1  -  'A2  -  1 ) ( Bt 2  -  Btl)/ICA2  -  1 )B, [  *  B( 2 

(  1 


I  ti  terms  of  A  and  rj ,  we  can  write  (31  a  nd  ( 9  ) . 
respectively,  as 


Pp/Pl 


1  -  (b[?  -  Bu) 


<  IBtl  *  B 1 2  -  -?(A-  -  1  )BnVB,  2  !/8*pj  (1 


and 


F2/Pl  “  *?  t  1(7  “  l)/yll(A/  -  l)Brj  + 


X  lir,?  -  l)(Bn2  +  Btl2)  -1  Bcl2  -  Bt?2|/8irBr2P1 


(1 


Subtracting  (13)  from  (14), 

1 7  _  1  +  (  ( 7  '  1 )/ f ) (  (A^  - 

X  U-)2  -  U(Bn2  :  Stl2) 

*  lBt?  -  +  Bt2 

/  8*p^  -  0 


we  obtain 

i )  Bt  i  +  Rti? 1 

1  Btl2  -  Bt2  ’  * 
-  2 (A2  -  l)Bn2 


/8fBt 


Dividing  equation  (IS)  bv  r)  -  1.  making  use  of 
(1?)  and  rearranging  the  result  as  a  cubic 
equation  in  B^*  we  ^ave 

V?’  +  D2Bt  2 2  f  DlBt?  f  D0  "  0  1  1 


The  coefficients  D-  are  fmvt ions  of  upstream 
rondi  t  i  oris 


cl' .  '  ;?  +  ;  TP/c  )  /  (  y 
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Alfven  Number  =  Un/(a  cos  9)  Un/(a  cos  0) 

Fig.  4.  For  each  given  value  of  the  jump  conditions  of  flow  properties 
such  as  the  total  pressure  ratio  P2*/Pj*  and  the  deflection  angle  6  for  slew- 
shocks  can  be  plotted  as  functions  of  A  and  0 ^  At  the  lower  limit  of  A  the 
deflection  angle  6-0  and  there  are  no  jumps  In  p,  T.  p,  B,  and  p*.  At  the 
upper  limit,  the  shocks  are  known  as  the  switch-off  shocks.  Slow  shocks  exist 
for  all  finite  values  of  fJ|.  When  0^  »  1,  slow  shocks  exist  in  a  verv  small 
domain  in  the  A,  # p - parame ter  space,  and  these  shocks  are  verv  weak. 


D0  -  -(7  ♦  1  ><A2  -  l)2Bn2Btl 
D,  -  (Sirypj  -  2Bn2  *  7B.,2  ♦  (7  -  1)(A2  - 
"  (A 2  -  1  > 

D,  -  ! 1  ♦  (2  -  7) ( A2  -  1 ) )B  • 

When  the  unsrream  conditions  arc  given,  wr  first 
.solve  (16)  for  B„9.  Next  ,  we  can  calculate  all 
other  jump  conditions.  The  solutions  with  the 
.shock  Alfven  number  A  greater  than  1  represent 
fast,  shocks,  and  the  solutions  with  A  less  than  1 
represent  slow  shocks.  The  cubic  equation  «'16) 
may  ha^e  one  or  three  real  roots,  however,  fot 
an v  given  upstream  conditions,  no  more  than  one 
root  produces  a  shock  solution  with  density  ratio 
r)  greater  than  1  . 

Algebraic  solut ions  of  MHD  shocks  havt  been 
ex'ensivelv  studied  bv  manv  authors  using  two 
general  methods  of  approach.  We  take  the  same 
approach  as  Lust.  [19V)j  using  only  the  upstream 
conditions  of  a  shock  as  the  input,  for  the 
calculation.  His  method  involves  the  solution  of 


a  cubic  equation  in  density  ratio.  The 
coefficients  of  the  cubic  equation  in  B^ 
obtained  in  (16)  are  much  simpler  than  their 
counterparts  obtained  by  Lust.  The  second  method 
involves  the  introduction  of  a  shock  parameter 
which  is  a  function  of  some  flow  conditions  on 
both  sides  of  the  shock  [Heifer,  1953.  Barer  and 
Ericson,  1959;  Kulikovskii  and  Lyubimov.  1^62; 
Lynn,  1966].  Insofar  as  the  downstream  conditions 
of  a  shock  are  completely  unknown,  as  in  the 
example  shown  in  the  next  section  for  the  slow  - 
fast  shock  system,  the  present  solution  is 
believed  to  be  a  more  direct  and  simple  method  of 
solution  for  MHD  shocks. 

As  shown  in  Figures  1  and  2,  for  each  given 
value  of  the  jump  conditions  can  he  plotted  as 

functions  of  A  and  9  ^ .  For  varving  values  of 
we  plot  the  constant  contours  of  the  total 
pressure  ratio  P2*''P\*  anci  the  deflection  angle  ft 
in  Figures  a  and  6  for  slow  shocks  and  fast 
shocks,  respectively.  The  existence  <>f  shock 
solutions  requires  that 

^nl  ^  l'sl 
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Alfven  N.jmber  -  Un/(a  cos  9)  AHven  Number  =  Un/(a  cos  0) 


Fig  5.  For  each  given  value  .»f  fi,  the  jura;  conditions  of  fast  shocks  ca:i 
plotted  as  functions  of  A  and  6  When  9  ^  -  0,  the  lover  limit  of  A  equals  1 
for  St  <  2/7  and  equals  7^/2  for  >  2/y.  Because  the  onset  ot  u„:  s;i  io:i 
from  slou  shock  to  fast  shock  can  take  place  onlv  when  the  lower  limit  of  A  : 
unity,  it  must,  occur  in  the  region  of  interplanet  trv  :.p.-ce  where  r  -  ! 


for  slow  shocks  and 

Unl  >  Cfl 

for  fast  shocks,  where  Cg  and  are  the  normal 
speed  (the  phase  speed)  of  the  slow  and  fast 
magnetoacoustic  wave 

Cs2  -  (a2  +  c2)/2 

-  I  (a2  +  c2)2  -  4a2c2cos2(? )  1/2/2  (17) 

and 

Cf2  -  (a2  +  c2)/2 

+  |(a2  +  cz)2  -  4a2c2cos2« ) 1/2/2  (18) 

where  c  is  the  gasdynamic  sound  speed,  c2  -  rp/r. 

From  (17)  we  obtain  that  for  given  0-^  and  6 j, 
slow  shocks  exist  in  the  interval 

1  +  yP^/2  -  {(1  +  7^;/2 ) 2  -  27/l1cos2  9  j  1 1  /2 

2  cos'^i 

<  A2  <  1 

These  limits  are  shown  in  Figure  6  for  various 
values  of  0-^.  At  the  lower  limit  the  deflection 
angle  <5-0,  and  there  is  no  jump  in  /> ,  T.  p.  3. 
or  p*.  At  the  upper  limit  the  deflection  angle  h 
equals  (i.e,,  "  0  and  ~  0),  and  the 


shocKS  are  known  as  ’’he  switch -off  shocks  When 
6  ]  -  0,  the  lower  limit  equals  yfi^/2  for  hi  ' 
and  equals  1  fo*  0  ^  >  2/y  When,  9  ^  -  ^0°,  the 
lower  limit  equals  (1  +  2 Note  ?haf  :  h.i  s 
limiting  value  of  A  is  less  than  1  for  all  f  * » 1 1  *■ 


values  of  0^.  Thus  slow  shocks  exis*  for  all 
finite  values  of  0^.  When  0^  >>  1,  slow  shocks 
exist  in  a  very  small  domain  in  the  A.0i* 


parameter  space,  and  these  shocks  are  vei 
At  the  lower  limit  of  A4  -  (1  +  2/y0^)~ 

«  90°,  the  group  speed  for  'he  slow  mod^ 
magnetoacoustic  wave  equals  (1/a"  +  l/c") 


v  weak 
when  f?  | 

-1/2 


which  is  the  magnetoacoustic  cusp  speed. 


From  (18)  we  find  that  for  given  0^  and  0., 
the  shock  Alfven  number  of  last  shocks  has  a 


lower  limit 


1  t  7/^/2  +  1(1  +  70;/2>2  -  2t^1cos2«1  1 1/2 

> - 

2  cos^| 

At  the  lower  limit  the  fast  shock  approaches  a 
fast  magnetoacoustic  wave,  the  deflection  angle  6 
equals  0,  and  there  is  no  jump  in  an.  flow 
property.  When  -  0,  the  lower  limit  equals  1 
for  0-^  <  2/y  and  equals  y  # ,  /  2  for  0^  >  2/y.  The 
onset  of  transition  from  slow  shock  to  fast  shock 
can  take  place  only  when  the  lower  limit  of  A  is 
unity  at  0  j  -  0 .  For  this  reason  the  onset  o  f 
transition  must  occur  in  the  region  of 
interplanetary  space  where  ft  <  2/y  nr  1.7. 

Let  us  now  examine  the  constant  contour  for 
the  deflection  angle  <S  in  Figure  N  for  -  0 
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Fig.  6.  The  configuration  of  a  slow-fast  shock 
system  predicted  by  this  study  for  the  transition 
of  traveling  forward  slow  shock  to  fast  shock. 

The  principal  components  of  the  system  are  a  slow 
shock  (S),  a  fast  shock  (F),  and  a  rotational 
discontinuity  (R) ,  as  shown  in  this  sketch.  The 
three  surfaces  of  discontinuity  intersect  along  a 
closed  loop  which  may  be  called  a  transition 
line.  The  strength  and  the  direction  of  each 
surface  of  discontinuity  make  sudden  changes 
across  the  transition  line. 


When  <  2/7,  the  shock  angle  jumps  from  0  to  a 
significantly  large  angle  for  A  slightly  greater 
than  1.  These  solutions  represent  switch-on 
shocks.  Switch-on  shocks  can  not  exist  in  the 
region  where  0  >  1.2. 

4,  The  Slow-Fast  Shock  System 

During  the  transition  the  slow-fast  shock 
system  has  a  complicated  configuration.  We  obtain 
a  numerical  solution  to  demonstrate  that  the 
principal  components  of  such  a  system  are  a  slow 
shock  (S),  a  fast  shock  (F),  and  a  rotational 
discontinuity  (R) ,  as  sketched  in  Figure  6.  The 
three  surfaces  of  discontinuity  intersect  along  a 
closed  curve  which  may  be  called  a  transition 
line.  The  strength  and  the  direction  of  each 
surface  of  discontinuity  make  sudden  changes 
across  the  transition  line.  Let  SI.  Rl,  and  FI 
denote  the  discontinuities  inside  the  closed  loop 
of  the  transition  line  and  S2,  R2,  and  F2  the 
discontinuities  outside  the  closed  loop.  At  the 
onset  of  the  transition  process  the  fast  shock 
and  the  rotational  discontinuity  are  very  weak 
discontinuities  As  the  system  moves  outward  from 
the  sun,  the  transition  line  moves  laterally 
across  the  field  lines  in  directions  facing  away 
from  the  noses  of  the  three  bow-shaped  surfaces 
of  d i scont inui ty ,  the  area  enclosed  by  the  closed 
loop  of  the  transition  line  expands  continuously, 
the  fast  shock  grows  stronger,  and  the  slow  shock 
becomes  weaker. 


The  proposed  slow-fast  shock  system  shown  in 
Figure  6  which  may  represent  a  physical  process 
for  the  transition  of  slew  shocks  to  fast  shocks 
in  the  inner  solar  wind  requires  that  the  flow 
properties  in  the  immediate  neighborhood  of  the 
transition  line  satisfy  all  jump  conditions  of 
MHD  discontinuity  surfaces.  The  method  of 
solution  is  to  treat  the  normal  directions  of  the 
four  shock  waves  and  the  rotational  angles  of  the 
two  rotational  discontinuities  a.c  six  unknown 
variables.  In  order  to  analyze  the  flow 
conditions  in  the  immediate  neighborhood  of  the 
transition  line  it  is  convenient  to  introduce  a 
coordinate  system  attached  to  the  transition  line 
T  as  shown  in  Figure  7.  We  choose  the  z  axis 
tangential  to  the  transition  line  and  the  z,x 
plane  parallel  to  the  magnetic  field  vector  B 
upstream  of  the  shock  system.  The  normal 
directions  of  all  surfaces  of  discontinuity  are 
parallel  to  the  x,y  plane. 

We  present  a  simple  numerical  example  in  which 
the  upstream  solar  wind  velocity  relative  to  the 
shock  system  U  is  parallel  to  the  x,y  plane.  The 
direction  of  the  x  axis  is  chosen  in  such  a  way 
that  Ux  is  positive  upstream  of  the  shock  system. 
To  have  the  shock  Alfven  number  for  the  slow 
shock  S2  less  than  1  and  that  for  the  fast  shock 
FI  greater  than  1,  Uy  must  be  positive.  The 
example  is  carried  out  here  using  the  following 
upstream  conditions: 


0  -  0.2 


U/a  -  1.2 

arctan(Uy/Ux)  -  30° 
arctan(UzA)x)  -  0° 

arctan(By/Bx)  -  0° 

arctan(Bz/Bx)  -  10° 


Let  us  first  identify  the  three  upstream 
parameters  for  the  slow  shock  S2:  0  is  given;  ti  < 
shock  Alfven  number  and  the  shock  angle  can  be 
calculated  once  the  normal  direction  of  S2  is 
given.  This  means  that  for  a  given  normal 
direction  of  S2  we  can  solve  the  cubic  equation 
(16)  for  B 1 2  and  calculate  the  magnetic  field 
vector  downstream  of  S2.  Then  we  use  (7)  to 
calculate  [Uc]  and  the  flow  velocity  vector 
downstream  of  S2.  B  and  Uz  downstream  of  S2  are 
not  zero  any  more.  The  next  step  is  to  calculate 
P2  and  P2  using  (11)  and  (13).  Therefore,  once 
the  normal  direction  of  S2  is  given,  we  can 
determine  all  flow  conditions  downstream  of  S2  by 
solving  a  system  of  nonlinear  algebraic 
equations.  Following  the  same  procedure,  once  the 
normal  direction  of  FI  is  given,  we  can  determine 
all  flow  conditions  downstream  of  FI. 

Since  the  shock  Alfven  number  of  Rl  and  R2 
must  be  unity,  this  condition  is  used  to 
determine  the  normal  direction  of  the  two 
rotational  discontinuities.  The  strength  of  a 
rotational  discontinuity  may  be  measured  by  the 
rotational  angle  of  the  tangential  magnetic  field 
vector  We  use  the  right-hand  rule  to  determine 
the  positive  direction  for  the  rotational  angle 
Now  we  can  calculate  all  flow  conditions 
downstream  of  Rl  and  R2,  respectively,  for  given 
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Fig. 7.  A  coordinate  system  attached  to  the 
transition  line  T  is  used  for  the  analysis  of  the 
flow  conditions  in  the  immediate  neighborhood  of 
T.  The  z  axis  is  tangential  to  the  transition 
line  and  the  z,x  plane  is  parallel  to  the 
magnetic  field  vector  B  upstream  of  the  shock 
system.  We  study  a  simple  case  in  which  the 
upstream  solar  wind  velocity  relative  to  the 
shock  system  U  is  parallel  to  the  x,y  plane.  The 
direction  of  the  x  axis  is  chosen  in  such  a  way 
chat  Ux  is  positive  upstream  of  the  shock  system. 
SI.  Rl,  and  FI  denote  the  discontinuities  inside 
Che  closed  loop  of  the  transition  line.  S2,  R2, 
and  F2  are  outside  the  closed  loop.  The  flow 
downstream  of  SI  and  F2  is  separated  by  a  contact 
discontinuity . 

rotational  angles.  At  this  point,  we  have 
calculated  the  flow  properties  U/a,  B,  p,  and  0 
downstream  of  Rl  and  R2. 

Once  the  normal  directions  of  F2  and  SI  are 
given,  we  can  calculate  the  three  upstream  shock 
parameters  (A,  0,  and  0)  for  each  shock.  Once 
again,  we  can  solve  the  nonlinear  algebraic 


system  (equations  (16),  (7),  (11),  and  (13))  to 
determine  the  flow  conditions  downstream  of  the 
two  shocks.  In  summary,  the  flow  conditions 
downstream  of  SI  and  F2  are  functions  of  six 
parameters:  the  normal  directions  of  the  four 
shock  waves  and  the  rotational  angles  of  the  two 
rotational  discontinuities.  The  six  parameters 
are  treated  as  six  unknown  variables  to  be 
determined  by  boundary  conditions  fot  the  flows 
downstream  of  SI  and  F2. 

The  flows  downstream  of  SI  and  F2  are 
separated  by  a  contact  discontinuity.  Across  the 
contact  surface  the  flow  must  satisfy  three 
conditions:  (1)  the  magnetic  field  B  be 
continuous,  (2)  the  thermal  pressure  p  be 
continuous,  and  (3)  the  velocity  vectors  be 
parallel.  The  three  conditions  are  actually 
represented  by  six  scalar  equations.  We  can 
numerically  solve  the  system  of  nonlinear 
equations  for  the  six  unknown  variables.  Note 
that  since  U  x  B  is  continuous  across  each  of  the 
MHD  shocks  and  rotational  discontinuities, 
conditions  1  and  3  guarantee  that  U  is  also 
continuous  across  the  contact  discontinuity. 
Following  the  procedure  outlined  above,  we  have 
obtained  a  numerical  solution  which  satisfies  the 
jump  conditions  across  the  four  shock  waves,  two 
rotational  discontinuities,  and  one  contact 
discontinuity.  Important  parameters  of  the 
solution  identified  with  each  shock  and 
rotational  discontinuity  obtained  from  this 
example  are  listed  in  Table  1. 

If  Bz  -  0  upstream  of  the  slow-fast  shock 
system,  Bz  -  0  everywhere  in  the  immediate 
neighborhood  of  T.  The  rotational  angles  are  zero 
for  both  Rl  and  R2.  The  problem  is  reduced  to 
solving  a  nonlinear  system  for  four  unknown 
directional  angles  of  the  shocks. 

The  analysis  of  the  flow  conditions  in  the 
immediate  neighborhood  of  the  transition  line 
suggests  that  during  the  transition  process  the 
large-scale  slow-fast  shock  system  can  have  a 
configuration  as  shown  in  Figure  6.  This  model 


TABLE  1.  An  Illustrative  Example  for  the  Slow-Fast  Shock  System 


Discontinuity  SI  S2  FI  F2  Rl  R2 


Normal 

Direction 

of  Each 

Discontinuity 

arc  tan(ny/nx) 

-48.8 

-44.0 

34.6 

37.2 

9.7 

-6.0 

Direction 

of  U  and  BBehind  Each  Discontinuity 

arc tan(U  /Ux) 

21.2 

32.3 

18.2 

21.2 

18.3 

31.6 

arctan(Uz/Ux) 

4.3 

1.1 

3.4 

4.3 

3.4 

1.5 

arctan(B  /Bx) 

-13.2 

-1.7 

-11 .4 

-13.2  ■ 

-11.4 

-2.9 

arctan(BzyBx) 

12.4 

9 . 7 

12.8 

12.4 

12.9 

1 0 . 2 

Rotation  Angle  of 

Bt 

-0.02 

7.06 

Jump 

Conditions 

of  the 

Four  Shocks 

*i 

Ai 

ei 

0.200 

0.200 

0.200 

0.322 

0.4  74 

0.467 

1.476 

1.663 

39.15 

44.89 

35.84 

41 . 19 

b2/b1 

0.976 

0.972 

1 .  196 

1  .  207 

&  1 

1  .  296 

1  .  283 

1.221 

1.239 

t2/ti 

1  . 195 

1  . 186 

1.173 

1  .  174 

P2/Pl 

1 . 548 

1.521 

1.432 

1  455 

h 

0.325 

0.322 

0.200 

0.322 

a2 

0.416 

0.412 

1 .  336 

1,493 

*2 

37.36 

4  3.19 

47.34 

5 1  , 44 
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for  the  transition  process  Infers  that  a 
traveling  forward  slow  shock  observed  In  the 
Inner  solar  wind  can  be  part  of  a  single  bow¬ 
shaped  slow  shock  surface  or  of  a  complicated 
slow-fast  shock  system  It  also  infers  that  the 
transition  process  provides  a  mechanism  which 
generates  rotational  discontinuities  in 
interplanetary  space.  The  presence  of  shock- 
rotational  discontinuity  sequences  in  the 
observational  data  (C.  F.  Kennel,  private 
communication,  1986)  possibly  supports  the 
existence  of  slow- fast  shock  systems  predicted  bv 
this  model. 
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Abs  t  ract ,  The  possible  existence  of  traveling 
forward  slow  shocks,  their  global  geometry  and 
their  transition  to  forward  fast  shocks  have  been 
discussed  in  a  recent  paper.  The  decrease  in  the 
Alfven  speed  at  increasing  heliocentric  distance 
causes  the  evolution  of  a  forward  slow  shock  into 
a  forward  fast  shock.  During  the  transition  the 
shock  system  consists  of  a  slow  shock,  a  fast 
shock,  and  a  rotational  discontinuity.  This  paper 
continues  to  discuss  three  aspects  about  the 
evolution  of  interplanetary  slow  shocks.  We  first 
present  a  survey  of  slow  shock  solutions  in  a 
three-dimensional  t\,8,0  parameter  space.  Here  A  - 
Un/(a  cos  8)  is  the  shock  Alfven  number,  Un  the 
normal  component  of  the  relative  shock  speed,  a 
the  Alfven  speed,  8  the  acute  angle  between  the 
shock  normal  and  the  magne: ic  field,  and  0  the 
ratio  of  the  thermal  pressure  p  to  the  magnetic 
pressure  B‘/8x.  In  a  region  where  the  plasma  0 
value  is  on  the  order  of  1  or  greater,  the  jumps 
in  thermodynamic  properties  across  a  slow  shock 
are  small  but  the  directional  changes  for  the 
magnetic  field  and  the  relative  flow  velocity  are 
not  necessarily  small.  On  the  other  hand,  in  the 
region  where  the  plasma  0  value  is  on  the  order 
of  0.1  or  less,  the  jumps  of  all  physical 
properties  across  a  slow  shock  may  vary  over  a 
wide  range  of  magnitudes.  Next,  we  discuss  that 
at  the  onset  of  the  transition  process  an 
interplanetary  slow  shock  smoothly  converts  to  a 
system  consisting  of  a  slow  shock,  a  very  weak 
rotational  discontinuity  (an  Alfven  wave;,  and  a 
very  weak  fast  shock  (a  fast  MHD  wave).  We  also 
show  that  during  the  transition,  the  amplitude  of 
fluctuations  in  flow  velocity  and  magnetic  fields 
are  large  in  the  turbulent  region  behind  the  fast 
shock . 

1.  Introduction 

Traveling  interplanetary  shocks  are  associated 
with  coronal  mass  ejections  and  flares.  When  the 
speed  of  the  ejected  material  relative  to  the 
ambient  solar  wind  exceeds  the  local 
magnetoacoustic  speed,  a  magnetohydrodynamic 
(MHD)  shock  front  forms  at  the  leading  edge  of 
the  compressed  ambient  plasma  shell.  The 
solutions  of  MHD  shocks  may  be  expressed  as 
functions  of  three  dimensionless  upstream 
conditions:  the  shock  Alfven  number  based  on  the 
relative  shock  speed  A  -  Un/(a  cos  8),  the  shock 
angle  8,  and  the  plasma  0  value.  Here  Un  is  the 
normal  component  of  the  relative  shock  speed,  a 
the  Alfven  speed,  $  the  acute  angle  between  the 
shock  normal  and  the  magnetic  field,  and  0  the 
ratio  of  the  thermal  pressure  p  to  the  magnetic 
pressure  B^/8jr.  Two  kinds  of  MHD  shocks,  the  fast 
and  the  slow  shocks,  may  exist.  A  shock  is  a  slow 
shock  if  the  shock  Alfven  number  A  is  less  than  1 
and  is  a  fast  shock  if  A  is  greater  than  1.  The 
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plasma  density,  the  thermal  pressure,  and  the 
total  pressure  (the  sum  of  the  thermal  and  the 
magnetic  pressures;  all  increase  across  an  MHD 
shock.  The  magnetic  field  and  the  shock  angle 
decrease  across  a  slow  shock  and  increase  across 
a  fast  shock.  Observat ions  of  slow  shocks  have 
been  reported  by  Chao  and  Olbert  [1970],  Burlaga 
and  Chao  [1971],  Richter  et  al.  [1985],  and 
Feldman  et  al.  [1985], 

In  a  recent  paper  [Whang,  1987]  (hereafter 
referred  to  as  paper  1),  we  studied  traveling 
interplanetary  forward  slow  shocks  and  their 
transition  to  forward  fast  shocks  in  the  inner 
solar  wind.  In  coronal  space,  the  Ait  en  speed  is 
of  the  order  of  1000  km/s;  the  shock  Alfven 
number  of  an  MHD  shock  in  the  inner  solar  wind 
can  be  less  than  1  and  the  shock  can  be  a  slow 
shock.  An  important  characteristic  of  slow  shocks 
(and  slow  MHD  waves)  is  that  they  are  upstream- 
inclined  [Sears  and  Resler,  1961],  The  large- 
scale  global  surface  of  a  forward  slow  shock 
should  have  a  bow-shaped  surface  with  its  nose 
facing  the  Sun.  Recently,  Hundhausen  et  al. 

[1987]  studied  the  outward  motion  for  coronal 
mass  ejections  observed  by  the  Skylab  coronagraph 
and  the  SMM  coronagraph.  Their  conclusion 
supports  the  possible  formation  of  forward  slow 
MHD  shocks  immediately  in  front  of  some  coronal 
mass  ejections.  They  pointed  out  that  in  a 
diverging  radial  background  field  a  large-scale 
traveling  forward  slow  shock  may  have  a  rather 
flat  configuration,  because  the  upstream  inclined 
nature  of  forward  slow  MHD  shocks  is  relative  to 
the  magnetic  field  geometry.  Thus,  the  observed 
flattening  of  the  tops  of  some  mass  ejection 
loops  support  the  possible  existence  of  traveling 
forward  slow  shocks. 

The  Alfven  speed  decreases  at  increasing 
heliocentric  distance  in  the  inner  solar  wind. 

The  variation  of  the  shock  speed  relative  to  the 
solar  wind  flow  is  not  well  understood.  Presume 
that  the  Alfven  speed  decreases  at  a  rate  faster 
than  the  relative  shock  speed  if  the  relative 
shock  speed  also  decreases  at  increasing 
heliocentric  distance.  Then  as  a  forward  slow 
shock  travels  outwards  from  the  Sun,  the  decrease 
in  the  Alfven  speed  causes  the  shock  Alfven 
number  of  a  forward  slow  shock  to  become  greater 
than  1  at  increasing  heliocentric  distances  and 
the  shock  evolves  from  a  forward  slow  shock  into 
a  forward  fast  shock.  Because  the  fast-mode  MHD 
waves  and  MHD  shocks  are  downstream- inclined,  the 
large-scale  global  surface  of  a  traveling  forward 
fast  shock  has  a  bow-shaped  surface  with  its  nose 
facing  away  from  the  Sun. 

Paper  1  suggested  a  possible  process  for  the 
transition  of  a  forward  slow  shock  to  a  forward 
fast  shock.  During  the  transition,  the  system 
consists  of  a  forward  slow  shock,  a  forward  fast 
shock,  and  a  forward  rotational  discontinuity. 

The  three  surfaces  of  discontinuity  intersect 
along  a  closed  transition  line.  At  the  onset  of 
the  transition  process  the  fast  shock  and  the 
rotational  discontinuity  are  very  weak.  As  the 
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Fig.  1.  This  figure  depicts  the  domain  in  which 
slow  shock  solutions  exist  in  the  three- 
dimensional  parameter  space.  The  domain  of  slow 
shock  solutions  in  the  A,0^  parameter  space 
shrinks  as  the  plasma  P-y  value  increases.  The 
major  change  takes  place  when  py  is  between  0.1 
and  X.  The  change  becomes  very  small  at  very 
large  or  very  small  py. 

system  moves  outward  from  the  Sun  the  transition 
line  moves  laterally  across  the  field  lines,  the 
area  enclosed  by  the  closed  loop  of  the 
transition  line  expands  continuously,  the  fast 
shock  grows  stronger,  and  the  slow  shock  becomes 
weaker.  Eventually,  the  forward  slow  shock 
diminishes  and  the  entire  system  evolves  into  a 
forward  fast  shock.  This  proposed  transition 
process  satisfies  the  necessary  condition  that 
the  flows  in  the  immediate  neighborhood  of  the 
transition  line  meet  all  MHD  jump  conditions.  On 
November  12,  1978,  ISEE  1  and  3  observed  a  high¬ 
speed,  quas i - paral le 1 ,  interplanetary  forward 
fast  shock  which  was  followed  by  a  rotational 
discontinuity  [Kennel  et  al.,  1984a,  bj.  This 
event  resembles  the  kind  of  transition  process 
described  here. 

In  this  paper,  we  will  discuss  three  aspects 
of  traveling  interplane tarv  forward  slow  shocks. 
We  will  first  present  a  survey  of  slow  shocks  in 
a  three  dimensional  parameter  space.  Then,  we 
will  discuss  the  onset  of  the  transition  process 
Finally,  we  will  discuss  the  turbulent  nature  of 
the  flow  field  behind  the  fast  shock  during  the 
transition  process. 

2.  A  Three-Dimensional  Parametric  Survey 
of  Slow  Shocks 


of  the  MHD  shock  equations  as  functions  of  three 
dimensionless  upstream  parameters:  the  shock 
Alfven  number  A,  the  shock  angle  By,  and  py.  The 
subscript  1  denotes  the  flow  condition  on  the 
front  side  (upstream)  of  a  shock,  and  2  denotes 
the  flow  condition  on  the  back  side  (downstream). 
The  plasma  density  p,  the  thermal  pressure  p,  and 
the  total  pressure  p*  (the  sum  of  the  thermal  and 
the  magnetic  pressures)  all  increase  across  an 
MHD  shock.  The  magnetic  field  and  the  shock  angle 
decrease  across  a  slow  shock.  The  decrease  in 
shock  angles  across  a  slow  shock  may  he  called 
the  deflection  angle  6.  In  this  section,  we 
present  the  solution  of  slow  shocks  in  a  three- 
dimensional  parameter  space,  the  A ,  -  space  as 

shown  in  Figures  1-3.  All  calculations  are 
carried  out  using  the  ratio  of  specific  heats  -y  - 
5/3  for  fully  ionized  plasma.  A  parametric  study 
of  MHD  slow  shocks  in  different  combination  of 
two  dimensional  parameter  space  has  also  been 
reported  by  Edmiston  and  Kennel  [1986’. 

Let.  us  first  examine  the  domain  in  which  slow 
shock  solutions  exist  in  the  three  dimensional 
parameter  space  (Figure  1).  Slow  shocks  exist  at. 
all  shock  angles  (from  0°  to  90°)  and  at  all  P 
values  in  a  certain  interval  of  the  shock  Alfven 
number  A, 

1  +  70j/2  -  ((1  +  7  x  /  2 ) 2  -  2v51<-'os2«1)1/2 
2  c  o  s  2  9 1 

<  A2  <  1 


DENSITY  RATIO 


A  =  Un/(a  cos  9) 


Fig.  2.  This  figure  presents  a  three-dimensional 
plot,  of  the  constant  contours  for  the  density 
ratios  across  MHD  slow  shocks.  The  jump  in 
density  across  a  slow  shock,  P^/^l’  var*es  over  a 
wide  range  of  magnitude,  from  verv  weak  to 
moderately  strong,  in  the  low  p  region  (where  py 
<  0.1).  On  the  other  hand,  onlv  weak 
discontinuity  in  density  across  a  slow  shock 
( p 2 / ^ i  <  1.5)  occurs  in  the  high  p  region  (where 
P}  '  >  1.0). 


In  paper  1,  we  obtained  a  direct  and  simple 
solution  for  the  jump  conditions  of  MHD  oblique 
shocks.  The  method  calculates  the  exact  solution 
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Fig.  3.  This  figure  presents  a  three-dimensional 
plot  of  the  constant  contours  for  the  pressure 
ratios  across  MUD  slow  shocks,  P2/pj/  We  can 
reach  the  same  conclusion  regarding  the 
discontinuity  in  pressure  as  functions  of  the 
plasma  0  value  from  Figures  2  and  3.  In  the 
coronal  space  where  the  0  value  of  the  solar  wind 
is  much  less  than  1,  the  change  in  thermodynamic 
properties  across  a  slow  shock  may  vary  over  e 
wide  range.  On  the  other  hand,  near  1  AU  where 
the  average  0  value  of  the  solar  wind  is  in  the 
order  of  1,  the  jump  in  thermodynamic  properties 
across  a  slow  shock  is  always  very  weak. 

The  lower  limit  is  equivalent  to  the  evolutionary 
condition  that 

Unl  >  Csl 

where  Cg  is  the  normal  speed  tthe  phase  speed)  of 
the  slow  magnetoacoustic  wave.  At  the  lower  limit 
the  deflection  angle  6-0,  and  there  is  no  jump 
in  p,  T,  p,  B,  or  p*.  When  9^  -  90°,  the  lower 
limit  equals  (1  +  2/i0^)~  .  When  -  0,  the 
lower  limit  equals  '10^/2  for  j. 3 <  2/7  and  equals 
1  for  0y  >  2/7.  The  slow  shocks  at  the  upper 
limit  of  A  -  1  are  the  switch-off  shocks  which 
turn  off  the  tangential  component  of  the  magnetic 
field  (i.e.,  $2  ~  0  and  S  -  0^).  The  domain  of 
slow  shock  solutions  in  the  A, Sj  parameter  space 
shrinks  as  the  plasma  0-^  value  increases  as  shown 
in  Figure  1  The  major  change  takes  place  when 
is  between  0.1  and  1.  The  change  becomes  very 
small  at  very  large  or  very  small  / 3^. 

Figures  2  and  3,  respectively,  present  the 
constant  contours  for  the  jumps  in  density  and 
pressure  across  slow  shocks,  namely  the  ratio 
P2^P\  or  P?/Pl'  as  functions  of  three  independent 
parameters,  A,  and  0^-  If  we  call  those  slow 
shocks  with  p^/P\  <  1.5  as  weak  discontinuities 
in  density  and  those  with  >  2.5  as  strong 


discontinuities  in  density,  then  from  Figure  2  we 
can  reach  a  conclusion  that  slow  shocks  covering 
a  wide  range  of  strength  for  the  discontinuity  in 
density  (from  very  weak  to  moderately  strong) 
-xist  in  the  low  0  region  (where  0-^  <  0.1).  On 
the  other  hand,  only  weak  discontinuities  in 
density  exist  in  the  high  0  region  (where  0-^  > 

1).  Next,  we  study  the  discontinuities  in 
pressure  for  slow  shocks.  From  a  plot  of  the 
pressure  ratios  in  Figure  3  we  can  reach  the  same 
conclusion.  In  the  coronal  space  where  the  0 
value  of  the  solar  wind  is  much  less  than  1,  the 
jump  in  thermodynamic  properties  across  a  slow 
shock  may  vary  over  a  wide  range.  On  the  other 
hand  near  1  At'  where  the  average  0  valu«  cf  the 
solar  wind  is  on  the  order  of  1,  the  jump  in 
thermodynamic  properties  across  a  slow  shock  is 
always  very  weak. 

The  jumps  in  other  physical  properties  offer 
different  ways  to  measure  the  effect  of  an  HMD 
slow  shock  on  the  plasma  and  fields.  The  constant 
contours  for  Ba/B^  at  0  -  0.2  were  plotted  in 
Figure  1  of  paper  1.  A  whole  range  of  changes  in 
Bj/Bi,  from  0  to  1,  may  take  place  in  the  domain 
of  slow  shock  solutions  at  0  -  0.2.  Across  an  MHD 
shock  the  magnetic  field  vector  abruptly  changes 
its  direction  by  the  deflection  angle.  The 
relative  flow  velocity  also  make  a  sudden  change 
in  its  direction  across  MHD  shock  by  an  angle 
approximately  in  the  similar  magnitude  as  the 
deflection  angle.  Constant  contours  for  the 
deflection  angle,  0  <  6  <  90°,  were  plotted  in 
Figure  9  of  paper  1.  All  field- related 
properties,  such  as  ^2^1  and  6,  change  over  the 
same  ranges  in  the  domain  of  slow  shock  solutions 
for  any  given  value  of  0 

The  plasma  0  value  is  on  the  order  of  0.1  in 
coronal  space,  it  increases  at  increasing 
heliocentric  distance.  In  the  region  where  the 
plasma  0  value  is  on  the  order  of  0.1  or  less, 
all  physical  properties  may  jump  across  a  slow 
shock  over  a  wide  range  magnitudes.  Near  1  AU 
where  the  plasma  0  value  is  on  the  order  of  1  or 
greater,  the  jumps  in  thermodynamic  properties 
across  a  slow  shock  must  be  small  but  the  jumps 
in  f ield- related  properties  are  not  necessarily 
small.  Traditionally,  we  measure  the  strength  of 
a  shock  by  the  jump  in  density  or  pressure  across 
the  shock  or  Po/Pj)-  Actually,  a  shock  may 

be  called  a  weak  shock  wave  only  when  the 
discontinuity  in  every  quantity  is  small.  A  slow 
shoex  with  small  discontinuity  in  tnermodynam ic 
properties  is  not  necessarily  a  weak  one.  When 
the  discontinuities  in  thermodynamic  properties 
are  small,  the  identification  of  a  slow  shock 
becomes  difficult  from  plasma  observations. 

3.  Onset  of  the  Transition  Process 

It  is  convenient  to  discuss  the  onset  of  the 
transition  process  using  an  AB  coordinate  system 
as  shown  in  Figure  A.  Here  A  is  the  shock  Alfven 
number  and  B  stands  for  the  ratio  of  the  magnetic 
field  B2/BJ.  We  may  call  the  region  in  which  both 
A  and  B2/B^  are  greater  than  0  and  less  than  1  a 
slow  region,  the  quadrant  in  which  both  A  and 
B2/B^  are  greater  than  1  the  fast  quadrant,  and 
the  point  at  which  both  A  and  B2/B^  are  unity  the 
Alfven  point.  The  domain  of  solutions  for  slow 
shocks  is  limited  to  the  slow  region  in  the  AB 
plane.  As  discussed  in  the  previous  section,  the 
domain  shrinks  as  the  plasma  0 |  value  increases. 
The  domain  of  solutions  for  the  fast  shocks  is 
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Alfven  Number  A  =  Un/(a  cos  8) 

Fig.  A.  The  upper  panel  plots  the  solution 
curves  for  MHD  shocks  with  0  >  2/7  and  the  lower 
panel  plots  the  solution  curves  with  0  <  2/7.  At 
the  onset  of  the  transition  process  the 
interplanetary  slow  shock  near  the  nose  smoothly 
converts  to  a  system  consisting  of  a  slow  shock, 
a  very  weak  rotational  discontinuity  (an  Alfven 
wave),  and  a  very  weak  fast  shock  (a  fast  MHD 
wave).  When  0  >  2/7,  fast  shock  solutions  do  not 
exist  near  the  Alfven  point  where  both  A  and 
B2/B^  are  unity.  Thus,  the  onset  of  transition 
must  occur  in  the  region  of  interplanetary  space 
where  0  <  2/7  or  1.2.  At  the  onset  of  transition 
the  slow  shock,  the  weak  rotational  discontinuity 
and  the  weak  fast  shock  move  at  the  same  speed 
equal  to  a  in  the  solar  wind  frame  of  reference. 

limited  to  the  fast  quadrant  of  the  AB  plane. 
Rotational  discontinuities  and  Alfven  waves  are 
represented  by  a  single  point,  the  Alfven  point 
o f  the  AB  p lane . 

A  forward  slow  shock  can  travel  in  the  inner 
solar  wind  having  a  bow-shaped  surface  with  its 
nose  facing  the  Sun.  A  large  pressure  jump 
bounded  by  the  discontinuity  surface  of  the 
forward  slow  shock  propagates  forward  in  the 
solar  wind  frame  of  reference.  On  the  curved 
surface  of  the  forward  slow  shock  the  shock 
Alfven  number  is  less  than  1  everywhere,  although 
It  varies  from  point  to  point.  The  slow  shock 
should  be  relatively  stronger  at  the  nose.  It 
becomes  weaker  on  the  flank  and  asymptotically 
approaches  slow  mode  MHD  waves.  Across  all  points 
of  this  single  surface  the  shock  solutions  are 
located  on  the  slow  region  of  the  AB  plane. 


As  the  forward  slow  shock  travels  outwards 
from  the  Sun,  the  decrease  in  the  Alfven  speed 
causes  the  shock  Alfven  number  to  reach  1  near 
the  nose  of  the  shock  surface,  the  transition 
from  a  forward  slow  shock  to  a  forward  fast,  shock 
begins  to  take  place.  At  the  onset  of  the 
transition  piuccss  the  interplanetary  forward 
slow  shock  smoothly  converts  to  a  system 
consisting  of  a  forward  slow  shock,  a  very  weak 
forward  rotational  discontinuity  (an  Alfven 
wave),  and  a  very  weak  forward  fast  shock  (a  fast 
MHD  wave).  On  the  AB  plane  the  mapping  of  the 
'low  sb  —  k  solutions  near  the  nose  lies  in  the 
immediate  neighborhood  of  the  Al^-'en  point  in  the 
slow  region.  As  shown  in  the  lower  panel  of 
Figure  A  the  shock  angle  approaches  0°.  The 
limiting  case  of  the  shock  is  a  gasdvnamic  normal 
shock  with  finite  strength.  The  slow  shock 
continues  to  carry  the  finite  jumps  in  presc  .... 
and  densitv  and  propagates  in  the  solar  wind 
frame.  The  rotational  discontinuity  maps  with  the 
Alfven  point.  On  the  AB  plane,  the  newly 
generated  weak  fast  shock  maps  with  points  in  the 
immediate  neighborhood  of  the  Alfven  point  in  the 
fast  quadrant. 

A  smooth  conversion  requires  that  at  the  onset 
of  transition  the  three  discontinuity  surfaces 
(the  slow  shock,  the  weak  rotational 
discontinuity  and  ’lie  fust  shock)  mow  at 

the  same  speed.  In  Figure  A,  the  upper  panel 
plots  the  solution  curves  for  MHD  shocks  on  the 
AB  plane  with  0  >  2/7  and  the  lower  panel  plots 
the  solution  curves  with  0  <  2/7.  When  0  >  2/7 , 
fast  shock  solutions  do  not  exist  near  the  Alfven 
point  in  Che  AB  p lane.  Thus,  a  smooth  conversion 
at  the  onset  of  transition  is  not  possible  when  0 
>  2/7  or  1.2.  The  onset  of  transition  must  occur 
in  the  region  of  interplanetary  space  where  0  < 
2/7  or  1.2. 

A.  Turbulent  Flow  behind  the  Fast  Shock 

Figure  5  shows  the  constant  contours  of 
deflection  angle  and  the  pressure  ratio  in  the 
domain  of  fast  shock  solutions  in  the  AS  plane  at 
a  given  0  value.  The  shock  Alfven  number  of  fast 
shocks  has  a  lower  limit 

1  +  10l/2  +  ((1  +  i0l/2)2  -  27,91cos2S1)1/2 

A2  > - - - — - 

2  cos2  9 j 

At  the  lower  limit  the  fast  shock  approaches  a 
fast  magnetoacoustic  wave,  the  deflection  angle  6 
equals  0,  and  there  is  no  jump  in  any  physical 
property . 

When  a  forward  transition  system  consisting  of 
a  slow  shock,  a  fast  shock,  and  a  rotational 
discontinuity  propagates  in  interplanetary  space 
with  <  1.2,  the  shock  Alfven  numbers  of  the 
fast  shock  near  the  nose  are  slightly  greater 
than  1  and  the  shock  angles  are  very  small.  From 
Figure  5  we  can  see  that  the  deflection  angles 
are  quite  significant  and  change  rapidly  for  fast 
shocks  in  this  neighborhood  of  the  parameter 
space.  The  shocks  signiricant.lv  change  the 
directions  of  the  magnetic  field  and  the  relative 
flow  velocity.  The  fast  shocks  in  this 
neighborhood  of  the  parameter  space  are  near 
switch-on  shocks.  Switch-on  shocks  and  near 
switch -on  shocks  may  exist  in  the  region  where 
<  1.2.  Small  random  fluctuations  in  plasma  flow 
and  magnetic  fields  on  the  front  side  of  the 
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Alfven  Number  A  =  Un/(a  cos  0) 


Fig.  5.  When  8  transition  system  propagates  in 
interplanetary  space  with  f)  <  1.2,  the  shock 
Alfven  numbers  of  the  fast  shock  near  the  nose 
are  slightly  greater  than  1  and  the  shock 
angles  are  very  small.  In  this  neighborhood  of 
the  parameter  space,  the  fast  shocks  signifi¬ 
cantly  change  the  directions  ot  the  magnetic 
field  and  the  relative  flow  velocity.  Small 
random  fluctuations  in  plasma  flow  and  magnetic 
fields  on  the  front  side  of  the  shock  are  con¬ 
siderably  amplified.  The  flow  becomes  very  tur¬ 
bulent  and  the  amplitude  of  fluctuations  in 
tiow  velocity  and  magnetic  fields  are  large  in 
the  turbulent  region  behind  the  fast  shock. 

shoe1/  are  considerably  amplified  across  the  near 
switch-on  shocks.  The  flow  becomes  very  turbulent 
behind  the  shocks.  Therefore,  during  the 
transition,  the  amplitude  of  fluctuations  in  flow 
velocity  and  magnetic  fields  are  large  in  the 
turbulent  region  behind  the  fast  shock.  When 
ISEE-1  and  -3  observed  a  possible  forward 
transition  system  consisting  of  a  forward  fast 
shock  ahead  of  a  rotational  discontinuity  on 
November  12,  1978,  an  intense  MHD  turbulence  was 
also  observed  in  the  region  between  the  two 
discontinuity  surfaces  [Kennel  et  al.,  1984b). 


5.  Summary 

In  this  paper,  we  discuss  three  aspects  of 
traveling  interplanetary  forward  slow  shocks:  (1) 
A  survey  of  slow  shocks  in  a  three-dimensional 
parameter  space,  (2)  the  onset  of  the  transition 
process,  and  (3)  the  turbulent  nature  of  the  flow 
field  behind  the  fast  shock  during  the  transition 
process . 

The  parametric  study  shows  that  in  the  region 
where  the  plasma  f)  value  is  of  the  order  of  0.1 
or  less,  the  jumps  of  all  physical  properties 
across  a  slow  shock  may  vary  over  a  wide  range  of 
magnitudes.  On  the  other  hand,  in  a  region  where 
the  plasma  f)  value  is  on  the  order  of  1  or 
greater,  the  jumps  in  thermodynamic  properties 
across  a  slow  shock  are  small  but  the  jumps  in 
field- related  properties  are  not  necessarily 
smal 1 . 

Next,  we  use  a  graphical  demonstration  to 
explain  that  the  onset  of  the  transition  process 
can  occur  in  the  region  of  interplanetary  space 
where  Q  <  2/y  or  1.2.  At  the  onset  of  the 
transition  process  the  interplanetary  slow  shock 
near  the  nose  of  the  slow  shock  smoothly  converts 
to  a  system  consisting  of  a  slow  shock,  a  very 
weak  rotational  discontinuity  (an  Alfven  wave), 
and  a  very  weak  fast  shock  (a  fast  MHD  wave).  The 


three  surfaces  of  discontinuity  move  at  nearly 
the  same  speed  (the  Alfven  speed)  in  the  solar 
wind  frame  of  reference. 

Finally,  we  show  that  during  the  transition, 
the  amplitude  of  fluctuations  in  flow  velocity 
and  magnetic  fields  are  large  in  the  turbulent 
legion  behind  the  fast  shock. 
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Abstract.  For  some  coronal  mass  ejections 
(CMEs),  their  interaction  with  the  ambient  solar 
wind  can  produce  a  forward- reverse  shock  pair. 

The  high-speed  mass  ejecta  compresses  the  plasma 
near  the  top  of  the  CME  on  both  sides  of  the 
tangential  discont inui tv  which  separates  the  CME 
plasma  from  the  ambient  solar  wind  plasma.  The 
front  of  the  compressed  CME  plasma  propagates  in 
the  reverse  direction  relative  to  the  ejecta 
flow,  it  may  steepen  to  from  a  reverse  slow 
shock  The  front  of  the  compressed  solar  wind 
plasma  also  propagates  in  the  forward  direction 
relative  to  the  ambient  solar  wind  and  it  may 
steepen  to  form  a  forward  shock.  The  forward- 

r cc,.i.'se  shock  [idu  _ _ _  .jltu  with  CMEs  moves 

outward  in  interplanetary  space  and  evolves  into 
a  pair  of  fast  shocks.  The  interplanetary 
manifestation  of  some  CMEs  is  pictured  as  a 
magnetic  cloud  accompanied  by  a  shock  pair:  a 
forward  shock  precedes  the  cloud  and  a  reverse 
shock  either  within  or  behind  the  cloud. 

1 .  Introduc  t ion 

The  possible  existence  of  forward  slow  shocks 
preceding  some  coronal  mass  ejections  (CMEs)  has 
been  suggested  by  Whang  [1987]  and  by  Hundhausen 
et  al.  [1987],  Whang  also  showed  that  transition 
from  a  traveling  forward  slow  shock  to  a  forward 
fast  shock  can  Cake  place  in  interplanetary 
space.  This  paper  suggests  that  a  reverse  slow 
shock  associated  with  CMEs  can  form  under  some 
circumstances,  and  it  shows  that  when  such  a 
shock  is  present  it  will  evolve  into  leverse  fast 
shock . 

The  solutions  of  magnetohydrodynamic  (MHD) 
shocks  may  be  expressed  as  functions  of  three 
dimensionless  upstream  parameters:  the  shock 
Alfven  number  based  on  the  relative  shock  speed  A 
-  Un/(a  cos  6),  the  shuck  angle  6,  and  the  plasma 
0  value.  Here  Un  is  the  normal  component  of  the 
relative  shock  speed,  a  the  Alfven  speed,  d  the 
acute  angle  between  the  shock  normal  and  the 
magnetic  field,  and  8  the  ratio  of  the  thermal 
pressure  p  to  the  magnetic  pressure  R^/8n.  Two 
kinds  of  MHD  shocks,  the  fast  and  the  slow 
shocks,  may  exist.  A  shock  is  a  slow  shock  if  the 
shock  Alfven  number  A  is  less  than  1,  a  fast 
shock  if  A  is  greater  than  1.  In  coronal  space, 
the  Alfven  speed  is  on  the  order  of  1000  km/s.  If 
a  forward  MHD  shock  forms  preceding  the 
compressed  ambient  plasma  shell  associated  with 
the  c Mr  shock  number  can  be  less  than 

1  and  the  shock  can  be  a  forward  slow  shock. 

The  plasma  density  p,  the  thermal  pressure  p, 
and  the  total  pressure  p*  (the  sum  of  the  thermal 
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and  the  magnetic  pressures)  all  increase  across 
an  MHD  shock.  The  magnetic  field  and  the  shock 
angle  decrease  across  a  slow  shock  and  increase 
across  a  fast  shock.  The  change  in  shock  angle 
across  an  MHD  shock  may  be  called  the  deflection 
angle.  Edmiston  and  Kennel  I  X 9 8 6  j  and  Whang 
f  1988]  have  studied  the  slow  shock  solutions  for 
the  jumps  of  physical  properties  as  functions  of 
various  combinations  of  dimensionless  parameters. 

The  plasma  0  value  plays  an  interesting  role 
in  the  jumps  of  thermodynamic  properties  across 
slow  shocks.  In  coronal  space  where  the  0  value 
of  the  solar  wind  is  on  the  order  of  0.1  or  less, 
the  discontinuity  in  thermodynamic  properties 
across  a  slow  shock  may  vary  over  a  wide  range  of 
magnitudes.  On  the  other  hand,  near  1  AU  where 
the  average  0  value  of  the  solar  wind  is  on  the 
order  of  1,  the  jump  in  thermodynamic  properties 
across  a  slow  shock  becomes  very  small.  The  jump 
in  f ield- related  properties  behaves  quite 
differently.  At  any  given  values  of  0,  in  the 
domain  of  solutions  for  slow  shocks  the  magnetic 
field  ratio  B0/B^  changes  over  a  whole  range, 
from  0  to  1,  and  the  deflection  angle  also 
changes  over  a  whore  range,  0  <  6  <  90°.  A  shock 
may  be  called  a  weak  shock  wave  only  when  the 
discontinuity  in  every  quantity  is  small.  A  slow 
shock  with  small  discontinuity  in  thermodynamic 
properties  is  not  necessarily  a  weak  shock. 
However,  when  the  discontinuities  in 
thermodynamic  properties  are  small,  the 
identification  o  .  a  slow  shock  from  plasma 
observations  becomes  difficult. 

2.  Forward-Reverse  Shock  Pairs 

We  would  like  to  investigate  the  dynamic 
interaction  between  the  coronal  mass  ejecta  and 
the  ambient  solar  wind.  Reference  should  be  made 
to  an  excellent  review  of  recent  work  on  CMEs  by 
Kahler  [1987],  The  association  between 
interplanetary  shocks  and  CMEs  has  been  clearly 
demonstrated  using  shocks  detected  at  Helios  1 
and  CMEs  observed  from  Solwind  coronagraph 
[Sheeley  et  al.,  1985],  With  very  few  exceptions, 
interplanetary  shocks  are  produced  by  major  CMEs 
and  are  confined  in  latitude  to  within  15°  of  the 
angular  extents  of  the  CMEs  measured  at  10  RQ 

The  number  of  slow  shocks  observed  in 
interplanetary  space  is  very  small.  Chao  and 
Olbert  [1970]  have  identified  two  forward  slow 
shocks.  Burlaga  and  Chao  [1971]  found  a  reverse 
a'd  a  forward  slow  shock,  Richter  ot  al.  [1981] 
reported  a  very  certain  identification  of  a 
forward  slow  shock. 

We  may  call  the  ejected  material  the  CME 
plasma  which  is  separated  from  the  solar  wind 
plasma  in  the  background  corona  by  a  tangential 
discontinuity.  The  large  momentum  of  the  high¬ 
speed  mass  ejecta  exerts  an  excessive  pressure  on 
the  ambient  plasma  shell  in  a  narrow  region  in 
front  of  the  CME.  This  causes  a  sudden  increase 
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Fig.  1.  The  direct  impact  of  high-speed  mass 
ejecta  on  the  ambient  solar  wind  compress  the 
plasma  near  the  top  of  CME  on  both  sides  of  the 
tangential  discontinuity.  This  interaction 
produces  a  forward  shock  preceding  the  CME  and  a 
reverse  shock  in  the  CME  plasma  closely  behind 
the  leading  edge. 


in  the  pressure  of  the  plasma  shell.  When  the 
velocity  at  the  front  of  the  CME  relative  to  the 
ambient  solar  wind  exceeds  the  local 
magnetoacoustic  speed,  a  forward  MHD  shock  may 
form  preceding  the  compressed  ambient  plasma 
shell.  The  nose  of  the  shock  should  be  at  some 
standoff  distance  from  the  top  of  the  CME.  The 
shock  is  a  forward  slow  shock  if  the  normal 
component  of  the  relative  shock  velocity  L'n  is 
greater  than  the  slow  mode  magne toacous t ic  speed 
but  less  than  a  cos  $.  The  shock  is  a  fast  shock 
if  U  ..  is  greater  than  both  the  fast-mode 
magne  toacous  t  i  c  speed  and  a  cos  6. 

Recent.lv.  Hundhausen  et  al.  [  1 9  8  7  J  studied 
CMEs  observed  from  the  coronagraphs  aboard  Skylab 
and  on  the  Solar  Maximum  Mission  satellite.  The 
observed  flattening  of  the  tops  of  some  mass 
ejection  loops  suggests  the  possible  existence  of 
traveling  forward  slow  shocks  immediately  in 
front  of  some  CMEs.  Their  study  of  the  observed 
speeds  for  the  outward  motion  of  CMEs  shows  that 
the  speed  at  which  a  CME  is  overtaking  the 
ambient  coronal  plasma  is  less  than  the  Alfven 
speed  but  greater  than  the  sound  speed  for  many 
(probably  most)  observed  CMEs.  This  also  supports 
the  possible  existence  of  forward  slow  shocks 

In  response  to  the  large  pressure  in  the 
compressed  ambient  plasma  shell,  the  CME  plasma 
immediately  behind  the  tangential  discontinuity 
is  also  decelerated,  compressed  and  heated.  This 
compressed  CME  plasma  is  confined  in  a  narrow 
leading  edge  region  bounded  by  the  tangential 
discontinuity  on  the  forward  side  and  by  a 
pressure  front  on  the  reverse  side.  The  pressure 
front,  propagates  in  the  reverse  direction 
relative  to  the  motion  of  the  ejected  material. 
Due  to  greater  wave  propagation  speed  on  the  side 
of  the  compressed  CME  plasma  at  the  pressure 
front,  the  pressure  profile  and  the  flow  speed 
profile  steepen  to  form  a  reverse  MHD  shock  as 
the  CME  moves  outward  from  the  Sun. 

The  ejected  CME  plasma  carrying  a  strong 
magnetic  field  is  a  low-0  plasma  in  which  the 
slow  mode  magnetoacoustic  speed  Cg  is  much  less 
than  the  Alfven  speed  a.  The  reverse  shock 
associated  with  CME  can  be  a  slow  shock  if  the 
normal  component  of  the  relative  shock  speed  Un 


is  greater  than  C  but  less  than  a  cos  *  Since  U 
is  the  relative  shock  speed  the  above  condition 
do e s  no t  i  m p  1  v  t ha t  t lie  ab s o  1  u  t e  w  1  oc  i  t  v  of  t  he 
ejecta  is  low.  We  do  not  vet  have  a  theory  to 
explain  whv  the  pressure  front  steepens  to  form  a 
slow  shock  rather  than  a  fast  shock  in  a  low-0 
plasma  in  the  coronal  environment.  The  possible 
observation  of  forward  slow  shocks  associated 
with  CME  by  Hundhausen  et  al.  ”1987]  provide;;  mu 
observational  support,  for  such  a  hvporh^sis.  fV. 
this  basis  we  estimate  that  the  reverse  shock  can 
be  a  slow  shock  whether  the  forward  one  is  a 
slow -mode  or  a  fast -mode  shock. 

Figure  1  shows  a  sketch  of  the  forward- reverse 
shock  pair  associated  with  CMEs  in  the  early 
phase  of  the  shock  pair.  A  large-  portion  of  t  he 
shock  surface  inside  the  ejecta  is  quasi 
perpendicular  shock  with  the  shuck  angle  * 
slightly  less  than  or  equal  90°.  Fdmiston  .ind 
Kenne  1  ;  1  9  0 '  ]  and  Whcng  '  i  v 8  / ,  198  8  sho  w*.  d  th.it 
slow  shocks  with  large  shock  angle  ?  exist  at  all 
0  values.  In  f‘hc  limit  of  d  -  ^)°,  the  deflect  i  or. 
angle  6  approaches  zero  and  t  he  ratio  of  t  h*. 
total  pressor?  p*->/?*i  approaches  unity  The 
thermal  pressure  alwavs  increases  and  the 
magnetic  pressure  alwavs  decreases  across  a  5;  low- 
shock.  but  the  two  effects  cancel  out  in  the 
limiting  case  of  ^  -  90°  Because  the  total 
pressure  p*  must  he  continuous  .*t  every  point 
across  a  tangential  discontinuity,  the  jump  in  p* 
across  the  reverse  slow  shock  must  be  maintained 
across  the  boundary  surface  of  the  CME  plasma  as 
the  shock  intersects  with  the  boundary.  This 
requires  that  the  reverse  shock  must  extend  to 
the  solar  wind  outside  of  the  tangential 
discontinuity.  The  dynamical  equilibrium  of  the 
solar  wind  flow  field  in  the  neighborhood  of  the 
ejecta  has  to  be  maintained  in  such  a  wav  that 
this  boundary  condition  be  satisfied.  In  Figure 
2.  we  use  an  r.t  diagram  to  demonstrate  the 
forward- reverse  shock  pair  near  the  fop  of  a  CME. 
As  mass  ejections  sweep  through  the  corona,  a 
tangential  discontinuity  separates  the  CME  plasma 
from  the  solar  wind  in  the  background  corona,  a 
reverse  slow  shock  propagates  in  the  region  of 
CME  plasma,  while  a  forward  shock  propagates  in 
the  solar  wind  ahead  of  the  CME.  The  dotted  lines 
indicate  the  paths  of  fluid  elements. 

Numerical  simulations  have  been  used  to 
demonstrate  the  possible  existence  of  forward 
shock  in  association  with  flares  bv  Hundhausen 
and  Gentry  '1969'  and  the  possible  existence  ot 
forward  reverse  MHD  shock  pairs  bv  Steinolfson  et 
al.  19  7  5;  and  Whang  [198  +  :.  The  evolution  of  a 
corotating  streams  also  leads  to  the  formation  of 
a  corotating  forward  -  reverse  shock  pair  [Dossier 
and  Fejer,  1963;  Sonne tt  and  Colburn.  1965;  Whang 
and  Chien.  1981].  The  region  bounded  by  a 
corotating  shock  pair  is  known  as  a  corotating 
interaction  region  (CIR).  Evolution  of  corotating 
shock  pairs  and  CIR  has  been  very  well  studied  in 
the  literature  (see  a  review  by  Burlaga  [1988] 
and  references  therein).  Co rotating  shock  pairs 
begin  to  form  at  near  1  AU  in  interplanetary 
space.  The  shocks  are  weak  near  1  AC  and  they 
become  quite  strong  in  the  outer  heliosphere.  The 
direction  of  a  corotating  stream  makes  a  spiral 
angle  with  the  heliocentric  radial  direction. 

This  effect  substantially  reduces  the  normal 
impact  of  the  momentum  of  a  high-speed  stream  on 
preceding  solar  wind  streams.  In  contrast,  the 
coronal  mass  ejecta  can  exert  a  diiect  impact  on 
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Fig.  2.  As  mass  ejections  sweep  through  t  he*  corona,  a  tangential  d :  scout  i  tv.:  i  t  v  s**par  a;  «-s 
the  CME  plasma  from  the  solar  wind  in  '.he  background  corona,  a  reverse  sl»*w  .shock  pmpa 
gates  in  the  region  ot  '-ME  plasma,  while  a  forward  shock  propagates  in  the  solar  wind 
ahead  of  the  CMK.  The  clotted  lines  indicate  the  paths  of  fluid  elements 


the  ambient  solar  wind  The  plasma  pressure  can 
build  up  very  rapidly  near  the  top  of  a  C ME  on 
both  sides  of  the  tangential  discontinuity  which 
separates  the  CM E  plasma  f i om  the  ambient  solar 
wind  plasma  Thus  the  torwarci  reverse  shock  pairs 
associated  with  CMEs  can  form  verv  rapidlv  in 
coronal  space. 

1.  Interplanetary  Evolution  of 
CME- Assoc iated  Shock  Pairs 

Two  important  evolutions  of  the  shock  pairs 
associated  with  CMEs  mav  occur  in  interplanetary 
spare  First,  a  CME  loop  mav  disconnect  from  the 
Sun  to  form  a  closed  magnetic  structure.  The 
disconnected  bubble  is  believed  to  manifest  as  a 
magnetic  cloud  in  interplanetary  space.  Second, 
slow  shocks  may  convert  to  last  shocks.  As  a 
result,  we  picture  the  interplanetary  consequence 
of  some  CMEs  as  a  magnetic  cloud  accompanied  by  a 
shock  pair:  a  Forward  fast  shock  precedes  the 
cloud  and  a  reverse  fast  si.cck  oirK^r  within  or 
behind  the  cloud. 

MacQueen  '1980]  studied  the  average  CME 
occurrence  rate  and  the  magnetic  flux  present  in 
each  CME  He  suggested  that  CMEs  do  not  carry  out 
a  distended  field  verv  far  to  interplanetary 
space,  but  rather  are  subject  to  a  reconnect  ion 
process  He  proposed  that  CME  loops  must 
magnetically  disconnect  from  the  Sun  and  form  a 
closed  magnetic  structure  The  disconnected 
bubble  continues  outward  into  interplanetary 
space 

Burlaga  et  al  1981]  analysed  the  magnetic 
field  cent igura t i on  behind  three  shocks.  Thev 
f  ound  a  system a t i c  con t i gu r a t i on  o f  the  m agne t i c 
field  and  called  if.  the  magnetic  cloud.  A 
magnetic  cloud  is  identified  as  an  interplanetary 
structure  in  which  the  magnetic  field  strength  is 
higher  than  average,  the  magnetic  field  direction 
rotates  monotonica 1 ly  through  a  large  angle,  the 
temperature  is  low,  and  the  plasma  $  is 


significantly  lower  than  1  Burlaga  and  B«.  h.inr.or 
19g2]  showed  that  clouds  persisted  to  the 
distances  of  2  to  9  Al"  and  estimated  that  the 
front  and  rear  of  a  cloud  expand  at  a  speed  ■{ 
the  order  of  half  the  Alfven  speed  Klein  and 
Burlaga  ]1982]  identified  *'o  clouds  in 
interplanetary  data  obtained  near  1  Al'  between 
19^7  and  19/8.  About  one  third  ot  them  are 
preceded  bv  a  shock,  and  these  clouds  appear  to 
he  moving  faster  than  the  ambient  solar  wind 
ahead  Thev  compared  the  physical  properties  of 
magnetic  clouds  with  those  of  CMEs  and  suggested 
a"  association  ul  some  clouds  with  disconnected 
magnetic  structures  of  CMEs  A  good  case  of  an 
association  between  a  CME  observed  bv  Sol  wind  and 
a  magnetic  cloud  observed  at  0.5*4  Al'  bv  Helios  1 
two  days  later  was  presented  bv  Burlaga  et  ai 
,1982]. 

The  Alfven  speed  decreases  at  increasing 
heliocentric  distance  in  the  inner  solar  wind  As 
a  forward  or  reverse  slow  shock  travels  outward 
from  the  Sun,  the  decrease  in  the  Alfven  speed 
causes  the  shock  Alfven  number  to  become  greater 
than  1  at  increasing  heliocentric  distances  and 
the  shock  evolves  from  a  slow  shock  into  a  fast 
shock.  The  transition  of  a  forward  slow  shock 
which  precedes  a  CME  to  a  forward  fast  shock  has 
been  discussed  bv  Whang  [1987]  The  reverse  slow 
shock  mav  also  evolve  into  a  reverse  fast  shock 
The  onset  for  the  transition  of  reverse  slow 
shocks  mav  take  place  at  the  tangential 
discontinuity  which  is  the  boundary  between  the 
CME  plasma  and  the  ambient  solar  wind,  inside  the 
boundary,  or  outside  the  boundary 

u ecu rrc n cy?  of  T r a ns_i  r  i on  I  ns  i d e  o i  Outside  o t 
CMEs 

I f  the  onset  of  transition  of  reverse  shocks 
occurs  either  inside  or  outside  of  the  boundary 
tangential  surface,  the  process  should  be  similar 
to  that  for  the  transition  of  forward  shocks 
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Transition  of  a  reverse  slow  shock  to  a  reverse  fast  shock 


Sun 


Fig  3  This  figure  depicts  the  process  for  the  transition  of  a  reverse  slow  shock  into 
a  reverse  fast  shock  The  normal  vectors  n  point  in  the  direction  of  plasma  tlow  relative 
to  the  shock  surface 

Figure  i  depicts  the  transition  of  a  reverse  slow  The  fast  shock  of  a  forward  or  reverse 

shock  to  a  reverse  fast  shock  The  vector  n  transition  svstem  is  a  near  switchon  shock, 

Heoor-oc  riiC  normal  cM  r«*rc  ion  pointing  in  the  across  which  small  r  .ndom  fluctuations  in  plasma 

direction  of  plasma  flow  relail”e  to  the  shock  tlow  and  magnetic  fields  on  the  front  side  of  the 

surface  The  solar  wind  plasma  enters  a  shock  shock  are  considerably  amplified  The  flow 

from  its  front  side  The  front  side  (the  upstream  becomes  verv  turbulent  behind  nearly  switchon 

side'  of  a  reverse  shock  is  the  side  closer  to  shocks  Therefore  during  the  transition,  the 

the  Sur  and  the  back  side  (the  downstream  side)  amplitude  of  fluctuations  in  flow  velocity  and 

is  the  side  further  awav  from  the  Sun.  The  magnetic  fields  are  large  in  the  turbulent  region 

surface  of  a  reverse  slow  shock  has  a  how-shaped  behind  the  fast  shock  Whang  1488 

surface  with  its  nose  facing  away  from  the  Sun  Kennel  et  al  ' 1984a,  V  reported  a  sequence 

and  the  surface  of  a  reverse  fast  shock  has  a  of  events  observed  from  ISEF.  1  and  1  that  seems 

bow-shaped  surface  with  its  nose  facing  the  Sun  to  resemble  the  forward  transition  svstem  On 

As  the  reverse  slow  shock  moves  outwards  from  November  1.’,  1478.  the  passage  over  the  lScr 

the  Sun,  the  decrease  in  the  Alfven  speed  causes  spacecraft  of  a  high  speed,  quasi -parallel . 

the  shock  Alfven  number  to  reach  1  near  the  nose  interplanetary  forward  fast  shock  with  l'  j  -  ?4() 

of  the  shock  surface,  the  transition  from  a  km  s'1,  9 ^  -  4l°,  a^  a  160  km  s'  and  3  "  -  1  .1  a 

reverse  slow  shock  to  a  reverse  fast  shock  begins  The  shock  is  followed  bv  a  strong  magnetic  field 

to  take  place  at  a  point  of  the  shock  surface  rotation  at  approx  i  mat  e  Iv  40  Rf.  downstream  of  the 

where  rhe  shock  angle  equals  0°  At  the  onset  of  shock,  and  an  extended  region  of  intense  MHD 

the  transition  process  the  reverse  slow  shock  turbulence  between  the  shock  and  the  magnetic 

smoothly  converts  to  a  system  consisting  of  a  field  rotation  Thev  reported  that  the  shock  was 

reverse  slow  shock,  a  very  weak  reverse  probablv  associated  with  the  pair  of  flares 

rotational  discontinuity  (an  Alfven  wave),  and  a  beginning  at  0048  and  0113  UT  on  November  10. 

very  weak  reverse  fast  shock  (a  fast  MHD  wave'  A  19"8 
smooth  conversion  requires  that  at  the  onset  of 

transition  the  three  discontinuity  surfaces  (the  Occurrence  of  Transition  at  the  Boundary  of  CMF.s 

slow  shock,  the  weak  rotational  discontinuity  and 

rhe  weak  fast  shock'  move  at  the  same  speed  Figure  1  shows  that  a  reverse  slow  shock  in 

whang  1987,  1988;  shows  that  rhe  onset  of  the  CME  plasma  extend  to  the  ambient  solar  wind 

transition  must  occur  in  the  region  of  outside  of  the  tangential  discontinuity  surface, 

interplanetary  space  where  8  '  i/i  or  i  /  it  rhe  relative  flow  velocity  across  the 

During  the  transition,  the  system  consists  of  tangential  discontinuity  is  less  than  the  Alfven 

a  reverse  slow  shock,  a  reverse  fast  shock,  and  a  speeds  on  both  sides,  the  flow  satisfies  the 

reverse  rotational  discontinuity.  The  three  criterion  of  Ke 1 vi n - He  1  mho  1 1 2  stability  The  flow 

surfaces  of  discontinuity  intersect  along  a  velocity  and  the  Alfven  speed  are  discontinuous 

closed  transition  line  As  the  system  moves  across  a  tangential  discontinuity  Rut  the  total 

outward  from  the  Sun.  the  transition  line  moves  pressure  p*  must  be  continuous  across  a 

laterally  across  the  field  lines,  the  area  tangential  discontinuity  and  the  magnetic  field 

enclosed  by  the  closed  loop  of  the  transition  vectors  and  the  relative  flow  velocities  on  both 

line  expands  continuously,  the  fast  shock  grows  sides  must  be  parallel  to  the  tangential 

stronger,  and  the  slow  shock  becomes  weaker  discontinuity  surface  As  the  shock  angle 

Eventually,  the  reverse  slow  shock  diminishes  and  decreases  across  the  reverse  slow  shock  inside 
the  entire  system  evolves  into  a  reverse  fast  the  CME,  the  boundary  surface  must  deflect  bv  a 

shock.  small  angle  at  the  location  where  the  reverse 
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slow  shock  intersects  the  tangential 
discontinuity.  The  flow  pattern  near  the 
deflection  point  of  the  boundary  in  the  ambient 
solar  wind  may  consists  of  a  reverse  slow  shock, 
a  reverse  rotational  discontinuity,  and  a  reverse 
fast  shock.  Depending  on  the  physical  properties 
outside  the  CME  boundary,  the  slow  (fast)  shock 
can  be  replaced  by  a  continuous  centered  slow 
(fast)  expansion  wave.  As  the  cluster  consisting 
of  the  cloud  and  the  shock  pair  moves  outward 
from  the  Sun,  the  decrease  in  the  Alfven  speed 
can  cause  a  possible  conversion  of  the  pattern 
for  shocks  and  discontinuities  near  the 
deflection  point  of  the  boundary  in  the  ambient 
solar  wind  from  the  pattern  of  a  reverse  slow 
followed  by  a  reverse  rotational  discontinuity  to 
that  of  a  reverse  fast  shock  preceded  by  a 
reverse  rotational  discontinuity 

4  Conclusions 

For  some  coronal  mass  eject  s  (CMEs).  the 
direct  impact  ot  the  high-speea  mass  ejecta  on 
the  ambient  solar  wind  can  substantially  compress 
the  plasma  near  the  leading  edge  of  the  CME  on 
both  sides  of  the  tangential  discontinuity  which 
separates  the  CME  plasma  from  the  ambient,  solar 
wind  plasma.  This  interaction  can  produce  a 
forward  shock  preceding  the  CME  and  a  reverse 
slow  shock  in  the  CME  plasma  closely  behind  the 
leading  edge.  As  the  forward- reverse  shock  pair 
moves  outward  in  interplanetary  space,  it  evolves 
into  a  pair  of  fast  shocks.  If  the  magnetic 
clouds  are  indeed  the  interplanetary 
manifestations  of  the  disconnected  magnetic 
structures  of  CMEs.  then  the  reverse  shock 
propagates  within  the  magnetic  cloud  as  the  cloud 
moves  outward  from  the  Sun  Eventually,  the 
reverse  shock  exits  the  cloud  to  became  detached 
from  the  magnetic  cloud.  Then  both  the  forward 
arid  the  reverse  shocks  propagate  in  the  ambient 
solar  wind  as  shown  in  Figure  4.  If  this  picture 
is  correct,  the  interplanetary  signature  of  some 
CMEs  at  large  heliocentric  distances  may  consist 
of  a  magnetic  cloud  wrapped  around  by  an 
interaction  region,  which  in  turn  is  bounded  by  a 
pair  of  fast  shocks  The  solar  wind  plasma 
continuously  adds  to  the  interaction  region 


Fig  4.  The  reverse  shock  propagates  within  the 
cloud  as  it  moves  outward  from  the  Sun. 
Eventually,  the  reverse  shock  exits  the  cloud  to 
become  detached  from  the  magnetic  cloud.  The 
cloud  is  now  wrapped  around  by  an  interaction 
region  which  is  itself  bound  bv  a  pair  of 
forward  reverse  fast  shocks. 


across  the  shocks  as  the  system  moves  away  I  rcrn 
the  Sun. 

Klein  and  Burlaga  [1982;  reported  that  about 
one  third  of  clouds  identified  near  1  AU  between 
1967  and  1978  are  preceded  by  a  shock  There  are 
two  examples  involving  the  observation  of  a 
reverse  shock  and  a  magnetic  cloud  R  P  Leppir.g 
et  al.  (private  communication,  1988)  have 
observed  one;  their  report  will  soon  be  ready  for 
journal  publication.  Another  event  involving  a 
reverse  shock  and  a  magnetic  cloud  was  observed 
on  August.  1982,  from  Voyager  7  Burlaga  et  a! 
1985],  The  observed  cloud  remained  stable  out  to 
10.3  AU  after  a  propagation  time  of  approximately 
a  month.  The  cloud  had  a  dimension  of  the  order 
of  1  AU  and  was  preceded  bv  two  shocks  and  an 
interface  and  followed  by  a  reverse  fast  shock  at 
about  0.1  AU  behind  the  cloud  This  observation 
probably  explains  that  for  a  long-lived  cloud, 
the  reverse  shock  had  enough  time  to  propagate 
through  the?  cloud,  exited  the  rear  of  the  cloud, 
and  appeared  closely  behind  the  cloud. 
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ABSTRACT 

For  some  coronal  mass  ejections  (CME~),  their  interaction  with  the  ambient  solar  wind  can 
0  produce  a  forward- reverse  shock  pair.  The  high-speed  mass  ejecta  compresses  the  plasma  near 

the  top  of  the  CME  on  both  sides  of  the  tangential  discontinuity  which  separates  the  CME 
plasma  from  the  ambient  solar  wind  plasma.  The  front  of  the  compressed  CME  plasma  propagates 
in  the  reverse  direction  relative  to  the  ejecta  flow,  it  may  steepen  to  from  a  reverse  slow 
shock.  The  front  of  the  compressed  solar  wind  plasma  also  propagates  in  the  forward  direction 
relative  to  the  ambient  solar  wind  and  it  may  steepen  to  form  a  forward  shock.  The  forward- 
reverse  shock  pair  associated  with  CMEs  moves  outward  in  interplanetary  space  and  evolves 
into  a  pair  of  fast  shocks.  The  interplanetary  manifestation  of  some  CMEs  is  pictured  as  a 
&  magnetic  cloud  accompanied  bv  a  shock  pair:  a  forward  shock  precedes  the  cloud  and  a  reverse 

shock  either  within  or  behind  the  cloud. 

INTRODUCTION 

i  he  solutions  of  MHD  shocks  /]/  may  be  expressed  as  functions  of  three  dimensionless  upstream 
conditions:  the  shock  Alfven  number  based  on  the  relative  shock  speed  A  -  U  /(a  cos  0).  the 
shock  angle  fi ,  and  the  plasma  0  value.  (Here  Un  is  the  normal  component  of  the  relative  shock 
4)  speed,  a  the  Alfven  .speed,  0  the  acute  angle  between  the  shock  normal^  and  the  magnetic  field. 

•  >nd  :i  the  ratio  of  the  thermal  pressure  p  to  the  magnetic  pressure  B' /8  ir.)  A  shock  is  a  slow 

shock  if  the  shock  Alfven  number  A  is  less  than  1  and  is  a  fast  shock  if  A  is  greater  than  1 
i  The  magnetic  field  and  the  shock  angle  decrease  across  a  slow  shock,  they  increase  across  a 

fast  shock. 


Hit-  association  between  interplanetary  .shocks  and  CMEs  has  been  cleariv  demonstrated  using 
Tanks  detected  at  Helios  1  and  CMEs  observed  from  Solwind  coronagraph  /2/.  In  1  ** 8 .  Whang 
first  suggested  the  possible  existence  of  forward  slow  shocks  preceding  some  CMEs  and  the 
transition  of  slow  .shocks  to  fast  shocks  /1,V-  The  theory  was  reported  at  the  XXVI  COS  PAR 
Meeting  at  Toulouse  on  July  ll)86.  Whang  explained  that  when  the  speed  of  ejected  CME  material 
relative  to  the  ambient  solar  wind  exceeds  the  local  magnetoacoustic  speed,  an  MHD  shock 
front  forms  at  the  leading  edge  of  the  compressed  ambient  plasma  shell  I  n  coronal  space,  the 
Alfven  speed  j  <;  of  the  order  of  1000  kin/s.  If  a  forward  MHD  shock  forms  preceding  the 
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fast  shock  of  a  transition  system  is  a  near  switch-on  shock  across  which  small  random 
fluctuations  in  plasma  flow  and  magnetic  fields  on  the  front  side  of  the  shock  are 
considerably  amplified.  The  flow  becomes  very  turbulent  behind  nearly  switch-on  shocks 
Therefore,  during  the  transition,  the  amplitude  of  fluctuations  in  flow  velocity  and  magm 
fields  are  large  in  the  turbulent  region  behind  the  fast  shock. 

Kennel  et  al.  /6,7/  reported  the  observation  of  a  rotational  discontinuity  preceded  bv  a 
forward  shock  from  ISEE  1  and  3.  On  12  November  1978,  the  passage  over  the  ISEE  spacecraft 
a  high  speed,  quasi  -  para  1 1  e  1 ,  interplanetary  forward  fast  shock  with  l'n]  -  240  km  s'^.  fj 
41°,  aj  =  160  km  s "  ^  and  =  1.14.  The  shock  was  probably  associated  with  a  pair  oi  flare 
which  occurred  some  48  hours  earlier.  The  shock  is  followed  b  ,■  a  strong  magnetic  field 
rotation  at  approximately  40  Rj.  downstream  of  the  shock,  and  an  extended  region  of  intense 
MHD  turbulence  hetween  the  shock  and  the  magnetic  field  rotation.  This  sequence  of  events 
observed  from  ISEE  1  and  3  seems  to  resemble  a  forward  transition  system. 


INTERPLANETARY  MANIFESTATION 


Two  important  evolutions  of  the  shock  pairs  associated  with  CMEs  may  occur  in  interplanetary 
space.  First,  a  CME  loop  may  disconnect  from  the  sun  to  form  a  closed  magnetic  structure.  Ti ,. 
disconnected  bubble  is  believed  to  manifest  as  a  magnetic  cloud  in  interplanetary  space. 
Second,  for  some  coronal  mass  ejections  (CMEs),  the  direct  impact  of  the  high-speed  mass 
ejecta  on  the  ambient  solar  wind  can  produce  a  forward  shock  preceding  the  CME  and  a  reverse 
slow  shock  in  the  CMF.  plasma  closely  behind  the  leading  edge.  As  the  forward  -  reverse  shock 
pair  moves  outwards  in  interplanetary  space,  they  evolve  into  a  pair  of  fast  shocks. 


MacQueen  /8/  suggested  that  CMEs  do  not  carry  out  distended  field  very  far  to  interplanetary 
space,  but  rather  are  subject  to  a  reconnection  process.  He  proposed  that  CME  loops  must 
magnetically  disconnect  from  the  sun  and  form  a  closed  magnetic  structure.  The  disconnected 
bubble  continues  outward  into  interplanetary  space. 


If  the  magnetic  clouds  are  indeed  the  Interplanetary  manifestations  of  the  disconnected 
magnetic  structures  of  CMEs  then  the  reverse  shock  propagates  within  the  magnetic  cloud  as 
che  cloud  moves  outward  from  the  sun.  Eventually,  the  reverse  shock  exits  the  cloud  to  became 
letached  from  the  magnetic  cloud  as  shown  in  Figure  3.  The  cloud  is  now  wrapped  around  bv  ar. 
.nteraction  region  which  is  itself  bound  by  a  pair  of  forward- reverse  fast  shocks. 


Burlaga  and  his  coworkers  /9-11/  Identified  a  magnetic  cloud  as  an  interplanetary  structure 
in  which  the  magnetic  field  strength  Is  higher  than  average,  the  magnetic  field  direction 
rotates  monotonically  through  a  large  angle,  the  temperature  Is  low,  and  the  plasma  S  is 
significantly  lower  than  1.  Clouds  may  persisted  to  the  heliocentric  distances  of  2  to  4  Al" 
They  also  suggested  an  association  of  some  clouds  with  disconnected  magnetic  structures  of 


Fig.  3  The  reverse  shock  propagates  within  the  c  1  oud  as  it  moves  outward  from  the  Sun. 
Eventually,  the  reverse  shock  exits  the  cloud  to  become  detached  from  the  magnetic  cloud.  Tlu 
cloud  is  now  wrapped  around  by  an  interaction  region  which  is  itself  hound  bv  a  pair  of 


CME  Associated  Forward-Reverse  Shock  Pairs 


(4)* 


forward- reverse  fast  shocks. 

CMEs.  Klein  and  Burlaga  /12/  reported  that  about  one  third  of  clouds  identified  near  1  AU 
between  1967  and  1978  were  preceded  by  a  shock. 

Two  recent  papers  reported  the  observations  of  CME  associated  forward- reverse  shock  pairs  at 
near  1  AU.  Lepping,  Ipavich  and  Burlaga  /13/  reported  that  a  large  magnetic  cloud  accompanie 
by  a  shock  pair  was  observed  on  October  31  and  November  1  of  1972  from  IMP  7.  The  plasma, 
magnetic  field  and  energetic  particle  data  observed  a  forward  shock  wave  followed  by  a  large 
tangential  discontinuity  and  a  reverse  shock.  The  tangential  discontinuity  was  the  forward 
boundary  of  the  magnetic  cloud.  The  reverse  shock  was  within  the  cloud.  The  flow  system  was 
associated  with  a  unusually  longlasting  solar  flare  which  occurred  some  49  hours  earlier. 
Gosling  et  al.  7H/  reported  that  they  have  identified  two  CME  associated  forward-reverse 
shock  pairs  at  1  AU  using  plasma  and  magnetic  field  data  from  ISEE  3  for  the  interval  from 
August  1978  through  February  1980.  The  shock  pairs  were  closely  associated  with  the 
propagation  of  transient  disturbances.  Strong  bidirectional  electron  heat  flux  events,  which 
they  identified  as  the  CME  driving  the  shocks,  were  detected  during  both  of  these 
disturbances.  In  the  stronger  of  the  two  disturbances  the  reverse  shock  was  found  within  tin- 
CME  and  was  separated  from  the  forward  shock  by  ~0.2  AU.  In  the  weaker  of  the  two 
disturbances  the  reverse  shock  had  propagated  entirely  through  the  CME  and  trailed  the 
forward  shock  by  0.3  -  0.4  AU. 

Burlaga  et  al.  / 1 5/  reported  another  event  observed  at  10.3  AU  from  Voyager  2  on  August  1987 
involving  forward,  reverse  shock  waves  and  a  magnetic  cloud  on  August,  1982  from  Voyager  2. 
The  observed  cloud  remained  stable  out  to  10.3  AU  after  a  propagation  time  of  approximately 
month.  The  cloud  had  a  dimension  of  the  order  of  1  AU  and  was  preceded  bv  two  shocks  and  an 
interface  and  followed  by  a  reverse  fast  shock  at  about  0.1  AU  behind  the  cloud.  This 
observation  indicates  that  for  a  long-lived  cloud,  the  reverse  shock  had  enough  time  to 
propagate  through  the  cloud,  exited  the  rear  of  the  cloud,  and  appeared  closely  behind  the 
cloud . 
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